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The interaction kernel in the Bethe-Salpeter equation for quark-antiquark bound states is derived 
from QCD in the general case where the quark and the antiquark can be of different flavors or the 
same flavor. The kernel is derived with the aid of the Bethe-Salpeter equations satisfied by the quark- 
antiquark four-point Green function. The latter equations are established based on the equations 
of motion obeyed by the quark and antiquark propagators, the quark-antiquark four-point Green 
function and some other kinds of Green functions in which the gluon field is involved. The interaction 
kernel derived is given a closed and explicit expression. Since the kernel is represented in terms of 

\f} • only a few-types of Green functions, it is not only convenient for perturbative calculations, but also 

f*"^ ' suitable for nonperturbative investigations. 
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I. INTRODUCTION 



o 

The Bethe-Salpeter (B-S) equation which was established early in Refs. [1, 2] is recognized as a rigorous approach 
to the relativistic bound state problem and has extensively been investigated in the period of more than half century 
[3-21]. The distinctive features of the equation are: (1) The equation is derived from quantum field theory and 
hence set up on the firm dynamical basis; (2) The interaction kernel in the equation contains all the interactions 
' taking place in the bound states and therefore the equation provides a possibility of exactly solving the problem 
of relativistic bound states; (3) The equation is elegantly formulated in a manifestly Lorentz-covariant form in the 
O^J ' Minkowski space which allows us to discuss the equation in any coordinate frame. However, there are tremendous 
difficulties in practical applications of the equation. One of the difficulties arises from the fact that the kernel in 
the equation was not given a closed form in the past. Usually, the kernel is defined as a sum of all B-S irreducible 
(or say, two-particle-irreducible) graphs. According to this definition, the kernel can only be calculated by the 
perturbation method. In practical applications, the perturbation series of the kernel has to be truncated in a ladder 
approximation. The ladder approximation has been proved to be successful in QED for studying the bound states 
Qh, formed by the electromagnetic interaction [22-26]. Nevertheless, it is not feasible in QCD for exploring multi-quark 
bound states because the confining force which must be taken into account in this case could not follow from a 
perturbative calculation, ft was remarked in Ref. [27] that " The approach using the Bethe-Salpeter equation has not 
led to a real breakthrough in our understanding of quark-quark force" . The reason is mainly due to that" we have 
no method for computing the kernel of the B-S equation" because we have not known the closed expression of the 
kernel which can be used to evaluate the confining force. In the previous application of the B-S equation to investigate 
the hadron structure, a phenomenological confining potential is necessarily introduced and added to the one-gluon 
exchange kernel so as to obtain reasonable theoretical results [27-29]. The confining potential was originally proposed 
from a nonperturbative computation, for example, the lattice gauge calculation for a special Wilson loop [30-33]. 

Since the B-S equation is an exact formalism for bound states and the B-S kernel contains all the interactions 
responsible for the formation of the bound states, the B-S kernel is supposed to be the most suitable starting point of 
examining the quark confinement. The question now is concentrated on whether a closed expression of the B-S kernel 
exists and can be derived? The aim of this paper is to give a definite answer to this question. It will be shown that 
opposite to the conventional concept that " The kernel K can not be given in a closed form expression" [34] , a closed 
expression of the B-S kernel for quark-antiquark bound states will perfectly be derived in this paper. The method 
of deriving the expression of the kernel is usage of various equations of motion satisfied by the qq four-point Green 
functions some of which involve the gluon field in them [19-21]. The kernel derived is compactly represented in terms 
of only a few types of Green functions. Therefore, the kernel is not only easily calculated by the perturbation method 
by employing the familiar perturbative expansions of the Green functions, but also suitable for nonperturbative 
investigations by the lattice approach by making use of the path-integral representations of the Green functions. 

The remainder of this paper is arranged as follows. In Sect. 2, we recapitulate the derivation of the B-S equation and 
show the definition of the B-S kernel. Section 3 serves to describe the procedure of deriving the B-S kernel. Section 4 
is used to derive the final expression of the kernel. In the last section, a summary and discussions arc presented. In 
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Appendix, the equations of motion satisfied by some Green functions which are necessary to use in the derivation of 
the B-S kernel are derived. 



II. DEFINITION OF THE B-S KERNEL 



To derive the B-S interaction kernel appearing in the B-S equation for qq bound states, we first show how the kernel 
is defined through Green functions. Before doing this, it is necessary to sketch the derivation of the B-S equation 
satisfied by the qq four-point Green function. For the qq bound states in which the quark and the antiquark may have 
same flavor, the four-point Green function is appropriately defined in the Heisenberg picture as follows [35] 



G(x 1 x 2 ;yi,y 2 ) 



rr 



(2.1) 



where ip(x) and ip c (x) are the quark and antiquark field operators respectively, ip(x) and ip (x) are their corresponding 
Dirac conjugates [5] 



V> c (x) = Cip (x), V C (x) = -?/> T (x)C- 



(2.2) 



here C = i-f 2 -f° is the charge conjugation operator, | } denote the physical vacuum states, T symbolizes the 
time-ordering product and N designates the normal product which is defined, for example, by 



N[^ a ( Xl )^(x 2 )} = T[ip a ( Xl )^(x2)} - (0+ |T[^ Q (x 1 )^(x 2 )]| 0-> . 



(2.3) 



It is emphasized here that the above normal product can only be viewed as a definition in the Heisenberg picture. 
With the definition shown in Eq. (2.3), the Green function in Eq. (2.1) can be represented as 



Q{xi,X2\yi,y2)aPpa = G(xiX 2 ; yi, yi)<xf)pa + Sp(X! - X 2 ) a f}S F {yi - y 2 ) f 

where 

G{x 1 x 2 ]y 1 ,y 2 ) al 3po- = (o + 
is the conventional qq four-point Green function, 



T{MxiW (x 2 )^ p ( yi )£(y 2 )} 



S£(xi - x 2 ) a p 



H 



T[^ a (xi)^(x 2 )] 



= S F {x\ - x 2 ) Q7 (C %p = S c F (x 2 - xi)p\Cx 



and 



in which 



and 



S* F (yi -V2)f 



0+ 



T{i/> p (yiW c M} 



= C crT S F (y 2 - y-Cjrp = (C ^pgSpiyi - y 2 ) SlJ 



S F (xi - x 2 ) 



(0+|T[^(x 1 )^ 7 (x 2 )]|0-) 



S c F (Vi - y2)5a = - (0+ T{^ 5 ( yi )£(y 2 )} 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2- 



(2.9) 



are the ordinary quark and antiquark propagators respectively [5]. It is clear that the propagators defined in Eqs. 
(2.6) and (2.7) are nonzero only for such a quark and an antiquark that they are of the same flavor. For the quark 
and antiquark of different flavors, the Green function defined in Eq. (2.1) is reduced to the conventional form shown 
in Eq. (2.5) since the second term on the right hand side (RHS) of Eq. (2.4) vanishes. In the case of the quark 
and antiquark of the same flavor, the normal product in Eq. (2.1) plays a role of excluding the contraction between 
the quark field and the antiquark one from the Green function. Physically, this avoids the qq annihilation to break 
stability of the bound state. It is pointed out that use of ?p c (x) other than ip{x) to represent the antiquark field in 
this paper has an advantage that the antiquark field would behave as a quark one in the B-S equation so that the 
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quark-antiquark equation formally is the same as the corresponding two-quark equation in the case that the quark 
and antiquark have different flavors. 

Now let us proceed to the derivation of the B-S equation for the qq bound system. This equation can be set up by 
operating on the Green's function Q{x\X2\ yi, 1/2) successively with the inverses of quark and antiquark propagators, 
i.e. the operator (id Xl —nii+ X) defined at the coordinate x\ and the operator (id X2 — m 2 + S c ) defined at the 
coordinate X2 where dx—j^d^, mi and are the quark and antiquark masses, and S and S c stand for the quark 
and antiquark proper self-energies which will be defined below. In Appendix, various equations of motion satisfied 
by several Green functions have been derived from the QCD generating functional. These equations of motion arc 
necessary for deriving the B-S equation. First, we show the equations of motion for the quark and antiquark two-point 
Green functions (the propagators). From Eq. (All) in the appendix A, we get the equation of motion satisfied by 
the quark propagator [5] 

[{id Xl -m 1 + E)5 F ] ap (a;i,yi) = <5 Qp ^ 4 (a;i - yi) (2.10) 

where 

(T,S F ) ap (xi,yi) = J d i zj:{x u zi) ai S F {z 1 -y 1 ) lp . . 

= (T»") Q7 Aj;(a;i|a;i,yi) 7/ , 

here 

(r a ") a7 = g(rT a ) ai (2.12) 
in which T a = A a /2 are the quark color matrices and 

A£(zi I zi,yi) 7 p = j (0+ |T[A°(a;i)^ 7 (a:i)^(yi)]| _ ) (2.13) 

in which A^(xi) are the gluon field operators. Similarly, from Eq. (B2) given in the appendix B, we can write the 
equation of motion for the antiquark propagator 

[(id X2 -m 2 + ^ c )S F W(x2,y2) = 5p a S 4 (x 2 - y 2 ) (2.14) 

where 

{H c S F ) p<T {x 2 , V 2) = J d 4 Z2^ c (x2, z 2 )pxS F (z2 - y 2 )xa 

= (T b ")e X Al»(x 2 \x 2 ,y 2 ) XrT (2 - 15j 



here 



t v = 91 »T b (2.16) 



with T b = — A b */2 being the antiquark color matrices and 

h cb {x2 \x2, V 2)xa = ^(0 + \T[A b ,(x 2 Wx(x2K(y2)]\0-). (2.17) 

Equations (2.11) and (2.15) just give the definition of the proper self-energies. The explicit representations of the 
self-energies can be found from the one-particle- irreducible decomposition of the Green functions A^(xi | xi,y±) and 

A^a* | aa.i/a) [5,36]. 

Here, we have a place to give general definitions of the Green functions involved in this paper for simplifying later 
statements. These Green's functions are 

A t'' b ^- d e( x ^ ■■■^ x f^y^ --^yi I xi,yi) ap , . 

= }(0+ | r[A-(x l )...A b (x,)A-(2 /fe )...A^(2 // )V Q (xi)^ p (yi)] | 0") 1 ' 8J 



^"vZJixii -,xf,yk, -,yi I x 2 , y 2 ) l3a _ 

= \(0+ I T[A-(x l )...Al(x 1 )AUy k ).--AUyiW (x2)£(y2)} | 0") 



(2.19) 



A ^«".e( a; *. -,Xj;yk, -,yi \ xi,x 2 ) af } 

- i(0+ | T[A-(x l )...At(x,)A-(y fc )...A^( 2/i )Va(x 1 )^(x 2 )] | 0") 



(2.20) 
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and 



K"- -vk'o <y - Xi ' ■ • ■ ' x 3 > v k > ■ • ■ ' f « I y 1 ' y -2 ) p° _ 

= i(0+ I T[A-(x l )...Al(x ] )A%(y k )...Al(y l )ij p (y 1 )£(y 2 )} | 0") 
G t'-' b ™' d e( x i' ■■■> x j'?yk, -,yi I x 1 ,x 2 ;yi,y 2 ) a p pa 

= (0+ | r[A»(^)-A fc ,(x J )A-(y fc )...A^(y ; )^(^i)^(^)V'p(yi)^(y2)] I o- 



■ e {x l ,...,x j ;y kl ...,yi \ x u x 2 ; yi, yi)oifi P c _ 
= (0+ | T{JV[AJ(a; < )...A*( a;j 0V'a(a:i)^(x2)]JV[^(i/i)<(|fi ! )AS(y fc )...Ag( w )]} I 0") 



(2.21) 



(2. 22) 



(2.23) 



where i,j.k, I = 1,2. The normal products in Eq. (2.23) are defined in the same way as that in Eq. (2.3). Hereafter, 
we no longer write individual representations of the Green functions encountered in later derivations . They can be 
read off from the above expressions. 

According to the derivations given in the appendices, we can write from Eqs. (A24) and (B4) the equations of 
motion obeyed by the Green function defined in Eq. (2.5) [5] 

[(id Xl -m 1 + ?,)G} a i3pa(xiX 2 ;yi,y 2 ) = <5 Qp <5 4 (a;i - yi)S F (x 2 - y 2 ) [ja 

+C a0 8 4 (xi -x 2 )S* F (yi -y 2 ) p<T - (T a ^) aj G a p (xi \ x u x 2 \y^y 2 ) 1&p( , (2.24) 
+ / d 4 z 1 T,(x 1 ,z 1 ) a7 G(z 1 ,x 2 ;y 1 ,y 2 ) 7 p pa 

and 

[(id X2 -m 2 + T, c )G] al3p a(xix 2 ;yi,y 2 ) = S /3a S 4 (x 2 - y 2 )S F (xi - yi) ap 

+C af! 5 4 (x 1 -x 2 )S* F (yi -y 2 ) pa - (T v )p\G h v (x 2 \ xi,x 2 ;yi,y 2 ) a \ pi7 (2.25) 
+ / d 4 z 2 Y;(x 2 ,z 2 ) l 3\G(x 1 ,z 2 ;y 1 ,y 2 ) a x p(7 

where the self-energy-related terms have been added on the both sides of the equations and G^(xi \ X\,x 2 ; yi,y 2 ) a p P cr 
was defined in Eq. (2.22) with one gluon field operator A^(xj) in it. It is noted that the second terms on the RHS 
of Eqs. (2.24) and (2.25) will disappear when the quark and the antiquark have different flavors. 

The equations of motion satisfied by the Green function defined in Eq. (2.1) can be written out by virtue of Eqs. 
(2.24) and (2.25). To do this, we need to introduce new Green functions which follows immediately from Eq. (2.23) 
as follows 

G*(xi | x 1 ,x 2 ;y 1 ,y 2 )aj3pa = G^(xi | x 1 ,x 2 ;y 1 ,y 2 ) a p pa 
+A* P a (x i | x 1 ,x 2 ) a ()S F {y 1 - y 2 ) P a 

where i = 1,2. On substituting the definitions in Eqs. (2.4) and (2.26) into Eqs. (2. 24) and (2.25), employing the 
equations in Eqs. (2.10) and (2.14) and the relations in Eqs. (2.6), (2.7), (2.11) and (2.15), it is not difficult to derive 
the following equations 

[(id Xl - mi + ^)Q] a p( r(xix 2 ; yi,y 2 ) = 5 ap 5 4 {xi - yi)S F (x 2 - y 2 )p a _ ^ 

-{^) ai Q p {xi | x 1 ,x 2 ;y 1 ,y 2 ) 1 p pa + J d 4 ziE(xi, z 1 ) ai Q{z 1 , x 2 ; y u y 2 )~ lf 3 pa : 

and 

[(id X2 -m 2 + T, c )g] af3p<y (xi x 2 ; yi,y 2 ) = S,3 a 5 4 (x 2 - y 2 )S F (x 1 - y{) ap 

xi,x 2 ;yi,y 2 ) aXpcr + J d 4 z 2 T, c (x 2l z 2 ) f 3\g(xi 1 z 2 ;yi,y 2 ) a x p(7 . 

In comparison of Eqs. (2.24) and (2.25) with Eqs. (2.27) and (2.28), we see that the second terms on the RHS of 
Eqs. (2.24) and (2.25), which are related to the annihilation between the quark and antiquark with the same flavor, 
are absent in Eqs. (2.27) and (2.28). That is to say, although the equations in Eqs. (2.27) and (2.28) suit to the both 
cases that the flavors of quark and antiquark are either different or the same, they formally are the same as those 
equations written in Eqs. (2.24) and (2.25) for the case that the quark and antiquark have different flavors. 

To derive the B-S equation, we need to operate on Eq. (2.28) with the operator( id Xl — mi+ E ) ( or cquivalently, 
operate on Eq. (2.27) with {id X2 -m 2 + S c )) 



[(id Xl - mi + S)(i9 X2 - m 2 + ^ c )Q] a i3 P a{xix 2 ]yi,y 2 ) 
= S l3a S 4 (x 2 - y 2 )[{id Xl - mi + T,)S F ] ap (xi - yi) 

-(r ")p\[{id Xl - mi + ?.)Gl\ a \ P o{x 2 | xi,x 2 ;yi,y 2 ) 
+ J d 4 z 2 Y, c (x 2 ,z 2 ) l3X [(id Xl -mi + ^)G] aXpa (xi,z 2 ; yi, 3/2)- 
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The first and third terms in the above can directly be evaluated by employing Eqs. (2.10) and (2.27). To compute 
the second term, we first write down the equation of motion for the Green function G b (x 2 \ Xi, x 2 ; J/i, y 2 ) which is 
given by supplementing the self- energy-related term to the both sides of Eq. (A25) 

[(id Xl - mi + ^)G b L/ } a \ pa (x 2 \ xi x 2 ; yi, y 2 ) = 8 ap S 4 (x 1 - ?/i)A£ b (x 2 \ x 2 , V2)\a 

+C aX S 4 (x 1 -x 2 )A* (x 2 I yi,y2) P a - (r a ^) Q7 G° b (xi,x 2 \ xi, x 2 ; yi, y 2 ) 7Ap(T ( 2 - 30 ) 

+ / d 4 ziS(a;i,zi) Q7 G^(x2 | z 1 ,x 2 ;yi,y2) 1 x P <7 

where A£ fc (x 2 \ x 2 ,y 2 )xcr, A*J'(x 2 | yi,y 2 )pa and G^ b (xi,x 2 \ xi, x 2 ; yi, y^Apcr are the Green functions which can be 
read off from Eqs. (2.19), (2.21) and (2.22) respectively. According to the relation shown in Eq. (2.26) and the 
following relation which follows from Eq. (2.23) 

Qf v {x\,X2 X\,X2\y\,y2)aJ)pa = G ab u (x 1 ,X 2 \ Xl, X 2 ; 2/1, y2)oc0pa , . 

+A™ b (x 1 ,X2\x 1 ,X2) a pS* F (y 1 -y 2 )p a 

where A*" fc (xi, X2 \ xi, x 2 )ai3 was defined in Eq. (2.20), and employing the equation of motion for the Green function 
A* b (x 2 I xi, X2) a p which is obtained from Eq. (A27) by introducing the self-cnergy-related term to the equation 



(2.32) 



i(id Xl - mi + £)A* b ] Q/3 (x 2 Xi, x 2 ) = -(r a ^) Q7 A*™(xi, x 2 |xi,x 2 ) 7/3 
+ j d 4 ziS(xi,zi) Q7 A* b (x 2 I z ll x 2 )~ / p 

it is easy to find from Eq. (2.30) an equation of motion for the Green function Gti x 2 | Xi, x 2 ; yi, y 2 ) such that 

[{id Xl - m± + Y,)Gl] a \po{x2 I xi,x 2 ;yi,y 2 ) = S ap 5 i (x 1 -y 1 )h cb {x 2 \ x 2l y 2 ) Xa 

+C aX S 4 (x 1 -x 2 )A* (x 2 I yi,y 2 ) pu - (r a ^) ai G p b u {x 1 ,x 2 \ x 1 ,x 2 ;y 1 ,y 2 ) 1 x prT (2.33) 
+ / d 4 z 1 Z(x 1 , z 1 ) ai g b (x 2 I zi,x 2 ;yi,y2)iX/><T- 

Upon inserting Eqs. (2.10), (2.27) and (2.33) into Eq. (2.29), we arrive at 

[(id Xl - mi + T,)(id X2 - m 2 + ^ c )G] a ppa{x\x 2 ; yi, y 2 ) ,^ ^ 

= S ap 8/3 a 6 4 (x 1 -y 1 )S 4 (x 2 -y2) + 'H 1 (x 1 ,x2;yi,y2) a i3pa 



where 



Hi{xi,x 2 ;yi,y 2 ) a/ 3pa = ^ Hi\xi, x 2 ; yi, y2) a p P a (2.35) 



in which 



H[ 1) (xi,x 2 - 1 yi, y2) a p P a = -(r aM ) Q7 J d A z2^ c {x2,z 2 )f}\Ql{x l \x l ,z 2 ;y 1 ,y 2 ) 1 x p(T , (2.36) 

(x 1 ,X2]yi,y2)aflpa = -(r^A / & Z\ S(xi , Zl ) al G b u (X2 I Z 1 ,X2]yi,y2) 7 Xpa, (2.37) 



Hi 3) (xi,x 2 ;yi,y 2 ) a(3(9(T = (r aAI ) a7 (r'" y ) ( gA^(xi,x 2 | xi, x 2 ; yi, t/2) 7 Ap<7, (2.38) 
Ti[ i) (x 1 ,x 2 ;yi,y2) a i3pa = <5 4 (xi - x 2 )(r'" / ) /3A GA Q A* fc (x 2 \ yi,yi) P a (2.39) 

W 1 5) (xi,x 2 ; y 1 , y 2 ) a p P a = J d 4 z 1 d 4 z 2 Y,{x ll z 1 ) ai Yf{x 2l z 2 ) 0X G{zi, z 2 ; yi,y 2 )jXpa- (2-40) 

From perturbative calculations or irreducible decompositions of the Green functions [1-5, 21], it may be found that 
the function Wi(xi, x 2 ; t/i, y 2 ) is B-S reducible and can be written in the form 



and 
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'Hi(xi,x 2 ;y 1 ,y 2 ) a0pr7 = J ' d 4 z 1 d 4 z 2 K(x 1 0jxQ(zi,z 2 ;yi,y2)jXpa (2.41) 

where K(x±, x 2 ; Zi, z 2 ) is precisely the B-S irreducible kernel. With the above expression for the function 
Tii(xi,x 2 ;yi,y 2 ), the equation in Eq. (2.34) can be expressed in a closed form 

[(id Xl - mi + T,)(id X2 - m 2 + ^Q^fjpaixi.x^y^y^ _ ^ 

= S ap SpaS 4 (xi -yi)S 4 (x 2 -y 2 ) + j d 4 z 1 d 4 z 2 K(x 1 ,x 2 ;z 1 ,z 2 ) af3lX g(z 1 ,z 2 ;y 1 ,y 2 ) 1 x pa . 

This just is the B-S equation satisfied by the Green function G{x\,x 2 ;yi,y 2 ). By making use of the Lehmann rep- 
resentation of the Green function G(xi,x 2 ;yi,y 2 ) or the well-known procedure proposed by Gell-Mann and Low [2], 
one may readily derive from Eq. (2.42) the B-S equation satisfied by B-S amplitudes describing the qq~ bound states 
[1-5] 

[(id Xl -mi + T,)(id X3 - m 2 + S c )xp ? ](a;i x 2 ) = j d i y x d A y 2 K{x u x 2 ; y u y 2 )xpdvi, Vz) (2-43) 

where 

Xp,{xix 2 ) = (0+ |JV[^(a:i)^ c (a? 2 )]| P?> (2.44) 

represents the B-S amplitude in which P denotes the total momentum of a qq~ bound state and <; marks the other 
quantum numbers of the state. The above equation can be written in the form of an integral equation when we 
operate on the both sides of the above equation with (id Xl — mi + T,)~ 1 (id X2 — m 2 + £ c ) _1 [l-5] 

Xp<;{xix 2 ) = j d 4 zid 4 z 2 d 4 yid 4 y 2 S F (xi - zi)S c F (x 2 - z 2 )K(z 1 , z 2 ; y u y^XPsiVu 2/2) (2.45) 



III. DERIVATION OF THE B-S KERNEL 



Beyond the perturbation method, the B-S kernel may be derived by starting from its definition shown in Eq. 
(2.41). One method of the derivation is usage of the technique of irreducible decomposition of Green functions. 
By this method, once the one and two-particle-irreducible decompositions of the Green functions contained in the 
function 7ii{xi, x 2 ; y\, y 2 ) are completed, it can be found that the function Hi(xi, x 2 ; yi, y 2 ) can really be represented 
in the form as written in Eq. (2.41) and, at the same time, the kernel will be explicitly worked out. Nevertheless, 
the expression of the kernel given in this way is much complicated ( We leave the detailed discussions in the future). 
Another method of deriving the kernel is to use B-S equations which describe variations of the Green functions 
involved in the function 7ii(xi,x 2 ;yi,y 2 ) with the coordinates y\ and y 2 . In this paper, we adopt the latter method 
and describe the procedure of the derivation in this section. Let us operate on the both sides of Eq. (2.41) from the 
right with the operator (i d yi + mi — d V2 + m 2 — S c ) such that 

/ d 4 z 1 d 4 z 2 K(x 1 ,x 2 ; Z!, z 2 ) a ^\[g(i d yi + mi - T,)(id y2 + m 2 - S c )] 7A(9CT (zi, z 2 ; yi, y 2 ) ^ ^ 
= Q(xi, x 2 ; yi,y 2 ) a /3 P a 

where 

Q{x 1 ,x 2 ;y 1 ,y 2 ) a p p(T = [Hi{id yi +m 1 - ^){id V2 +m 2 - T 1 c )] afjpiy (x l , x 2 ; y u y 2 ). (3.2) 

As seen from Eqs. (3.1) and (3.2), to derive the B-S kernel, it is necessary to use the B-S equations with re- 
spect to yi and y 2 for the Green function Q{x\, x 2 ; y\, y 2 ) and the other Green functions appearing in the function 
Hi(x\,x 2 ; yi,y 2 ). Analogous to the procedure stated in the preceding section, the above-mentioned B-S equations can 
be derived with the aid of equations of motion with respect to yi and y 2 for those Green functions. First, we show 
the equations of motion obeyed by the quark and antiquark propagators 

[S F (id yi +mi - £)] Qp (xi,yi) = -<5 Q(9 <5 4 (xi - y x ) (3.3) 

and 

[S F (id V2 +m 2 - E c )] /3ct (x 2 , y 2 ) = -5^5 4 (x2 - 2/2) (3.4) 
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which are given by Eqs. (C2) and (D2) respectively with the source J being set to vanish. We note here that the 
self-energies, as easily proved, are the same as those defined in Eqs. (2.11) and (2.15). 

Next, we derive the B-S equation satisfied by the Green function G(xi,x 2 ;yi,y 2 ). For this purpose, we first write 
down the equation of motion describing variation of the conventional Green function G{x\,x 2 \yi,y 2 ) with y\ which 
is obtained from Eq. (C5) by setting J = and adding the term related to the self-energy 

[G(id yi + m 1 - T,)] a/3pa (x 1 ,X2;yi,y2) = S ap S 4 (x 1 - yi)S F (x 2 - y 2 )fja 

-Cp (r 6 4 (yi-y 2 )Sp(x 1 -X2) a + G^(y 1 | x lt x 2 : yi, y 2 ) a p T a(r afJ -) T p (3.5) 
- / d 4 z 1 G(x 1 ,x 2 ; z 1 ,y 2 ) al3Ta Y,(z 1 ,y 1 ) Tp 

where G p (yi \ X\,x 2 \y\,y 2 ) is the Green function with a gluon field at y\ in it whose expression can be read from 
Eq. (2.22). Substituting in Eq. (3.5) the relation given in Eq. (2.4) and the following relation which follows from Eq. 
(2.23) 

GpiVi I x x ,x 2 \y x ,y 2 ) aPpa = G°(yi \ x 1 ,x 2 ;y 1 ,y 2 ) aPprT + Sp(x! - x 2 ) af jA* p a ( Vl | yi,y 2 ) p<7 (3.6) 
and employing the equation 

{S* F (id Vl + mi - S)U(»i,|ft) = C p J 4 ( yi - y 2 ) (3.7) 
which is given by Eqs. (C8)-(C10), it can be found 



[G(id yi + mi - T,)] a/3prT (xi,x 2 ;yi,y 2 ) = -5 ap 5 4 (xi - yi)S F (x 2 - y 2 )p a 
+%{y\ I x 1 ,x 2 ;y l7 y 2 ) al3Ta (T atl ) Tp - J d A ziG(xi,x 2 ; z l7 y 2 ) a , 3Ta T l (z 1 ,y 1 ) ap . 



(3.8) 



In order to derive the B-S equation, we need to operate on the both sides of the above equation with the operator 

{i*d y2 +m 2 - S c ), 

[G(i*d Vl +mi - Y)(id V2 +m 2 - Z c )] a /3 P a(xi,x 2 ; yi,y 2 ) = 5 ap 5 fja 8 i {x l - yi)S i (x 2 - y 2 ) 

+ [G P l(^d y2 +m 2 - £ c )W(yi | x 1 ,x 2 -y 1 ,y 2 )(T a ^) Tp (3.9) 

- / d 4 z 1 [G(id V2 +m 2 - S c )] Q(3TCT (a;i,a;2;2;i,y2)S(zi,2;i) T(9 

where Eq. (3.4) has been used. For computing the third term in the equation, we start from the equation of motion 
which is obtained from Eq. (D8) with setting J — and supplementing the self-energy-related term 

[G(id y2 +m 2 - £ c )] Q/3 p CT (xi,2;2;2/i,y2) = S PrJ S 4 (x 2 -y 2 )S F {xi -y\) ap 

-C pa 5 4 {yi -y 2 )S F (xi - x 2 ) af3 + G b v (y 2 | xi,x 2 ;yi,y 2 ) af 3 p s (FV (3.10) 
- / d 4 z 2 G(x 1 ,x 2 ; yi,z 2 ) al3p s^ c (z 2 ,y 2 ) Sa . 

By applying the relations shown in Eqs. (2.4) and (3.6) to Eq. (3.10) and considering the following equation 

\S F (id y2 +m 2 - T, c )] pa {y llV2 ) - C pa S 4 ( yi - y 2 ) (3.11) 

which is given by Eqs. (D14)-(D16) and the relations represented in Eq. (2.15), one can get such an equation that 

[G{id V2 +m 2 - Z c )} afjpry (xi,x 2 ; y 1 , y 2 ) = -5p a 5 4 {x 2 - y 2 )S F (xi - yi) ap 
+Gt(V2 I x 1 ,x 2 ;y 1 ,y 2 ) a(3p5 (T ) Sa - J d 4 z 2 G{x!,x 2 ; yi, z 2 ) a p pS T, c (z 2 ,y 2 ) Scr . 

To calculate the second term in Eq. (3.9), let us operate on the both sides of Eq. (3.6) with i = 1 by the operator 

{idy 2 +m 2 - E c ) 

[G p (id y2 +m 2 - E c )] Q/3(9a (yi | x 1 ,x 2 ;y 1 ,y 2 ) 

= [G%(i 9 V2 + m 2 - ^ c J\ a 0pa(yi \x- l ,x 2 ;y 1 ,y 2 ) (3.13) 
+5^(a;i - x 2 ) al3 [A* p a (id V2 +m 2 - S c )] pCT (j/i | yi,y 2 )- 

The first term can be calculated by virtue of the equation of motion derived in Eq. (D9). With addition of the 
self-energy-relevant term to Eq. (D9), we have 
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[G^(id V2 + m 2 - S c )] ct(3(9(T (2/i | xi, x 2 ; 2/1, 2/2) = -5 p<y 5 i {x 2 - y 2 )A p [(y 1 | xi,yi) c 



-C p(J <5 4 (yi - y 2 )A* a (yi | ari, a; 2 ) a/ 3 + Gf v {y u yi \x 1 ,x 2 - 1 y 1 ,y 2 ) a(ip i,{Y hv ) ia (3-14) 
-/ d 4 z 2 G°(j/i I xi,x 2 ;yi, z 2 ) a fj p sT, c (z 2 ,y 2 )s <J 

where the Green functions have been defined in Eqs. (2.18), (2.20) and (2.22). In accordance with Eq.(2.23), one can 
write 

Qf v {y\-,yi I x 1 ,x 2 ;yuy2)a0 P a- = Gf v {yi,y 2 \x\,x 2 \y\,y 2 ) a&pa 
+S* F (x 1 - x 2 ) af3 All (y 1 ,y 2 \y 1 ,y 2 ) pa 

where A*" fc (j/i,y 2 | yi,y2) P a was defined in Eq. (2.21). When the relations given in Eqs. (3.6) and (3.15) are 
substituted into Eq. (3.14), we are led to 

[G°(i*dy 2 +m 2 - S c )] ct( g p(T (yi | x\,x 2 ;y\,y 2 ) = -5p a 5 4 (x 2 - 2/2) A£ (2/1 | xi,yi) ap 
-C pa 5 4 (y 1 - y 2 )A* p a {y 1 | x\, x 2 ) a0 + Gf v {yi, 2/2 \ x 1 ,x 2 ;yi,y 2 ) a p p s(T V ) S(7 
- J* d 4 22^(2/1 I a;i,x 2 ;yi,z 2 ) Q)3 p < 5i; c (z 2 ,2/ 2 ) ( 5 (T - 5>(xi - x 2 ) Q/3 A*" (2/1,2/2 | 2/1, y2) P s(T 
+S* F (x 1 -x 2 ) af} J d 4 z 2 A* a (yi I yi,z 2 ) pS T l c (z 2 ,y 2 )sa- 

The second term in Eq. (3.13) is directly determined by the following equation of motion 

2/2 + m 2 ^ C )]pa(yi I 2/1)2/2) — A^^ (j/1,2/2 I 2/1) 2/2) p5 (r )5<r ^ -^y-j 

-/ d 4 z 2 A p a (y 1 I 2/i,z 2 ) p ,5E c (z 2 ,2/ 2 ) ( 5 (J 

which is given in Eq. (D17) with setting J = and adding the term associated with the self-energy. On inserting 
Eqs. (3.16) and (3.17) into Eq. (3.13), it is easy to find 

[G p (i < d y2 + m 2 - S c )] Q/3pCT (2/i I x 1 ,x 2 ;y 1 ,y 2 ) = -8 0tT 5 4 (x 2 - y 2 )A p (yi \ xi,yi) ap 
C pa S 4 (j/i - y 2 ) A* a (yi I Xi,x 2 ) a/ 3 + £^(2/1,2/2 l^i) ^2; 2/1) 2/ 2 W(r% CT (3-18) 



(3.16) 



-/ d 4 Z20 p (2/i I x 1 ,x 2 ;y 1 ,z 2 ) a p p sZ c {z 2 ,y 2 )s<,. 
Substitution of Eqs. (3.12) and (3.18) in Eq. (3.9) yields 



[G(id Vl +mi - T,)(id V2 +m 2 - S c )] Q/3p(T (a;i, x 2 ; y 1: y 2 ) ^ ig x 

= 8 ap 5p a 5 A {xi -yi)5 4 {x 2 -y 2 ) + H 2 {x\,x 2 ;yi,y 2 ) a p pa 



where 



5 

H 2 {x 1 ,x 2 ;y l ,y 2 ) OLl j pa = ^U 2 ] {xi,x 2 ;yi,y2) a f3pa (3.20) 



in which 



H { 2 \x 1 ,x 2 ;y 1 ,y 2 ) a0p<7 = - J d 4 z 2 G p {yi \ x x , x 2 ; 2/1, z 2 ) al3T sT, c (z 2 , 2/2)<5<x(r QAI ) rp , (3.21) 
H^ 2) (xi,a; 2 ;2/i,2/2)a/3 P <T = - / d 4 z 1 g b v (y 2 \ x ll x 2 \z ll y 2 ) al 3 T sYs{z 1 ,y 1 ) Tp (T hu )s IJl (3.22) 



H 2 ) (x l7 x 2 ;y 1 ,y 2 ) a f ipa = G p lb „(y 1 ,y 2 | Zi, x 2 ; 2/1, 2/ 2 )a/3T<5(r aA ') Tp (r >y ) ( 5 (T , (3.23) 
Hftxi,^;!/!,^)^ = <5 4 (2/i -y2)A; a ( yi \ x u x 2 ) a p{CY a ») ap (3.24) 



and 
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H^\x 1 ,x 2 ;yi,y2)a0p<7 = J d i z 1 d 4: z 2 g(x 1 ,x 2 ; Zi, z 2 ) a/3T s^{z 1} y 1 ) Tp T, c (z 2 ,y 2 )sa- (3.25) 
By applying the B-S equation given in Eq. (3.19), Eq. (3.1) can be written as 

K(xi,X2;yi,y2)af3pa = Q(xi,X2;yi,y2)a0pa (3 26) 

- J d i zid i z 2 K{xi, x 2 ; zi, z 2 ) a p T sH 2 {zi, z 2 ; yi, y 2 ) T s pa - 

To reach a closed expression of the B-S kernel, it is necessary to eliminate the kernel appearing in the second term 
on the RHS of Eq. (3.26). For this purpose, we operate on the both sides of Eq. (2.41) with the inverse of the Green 
function G{x\,x 2 ;yi,y 2 ). With this operation , noticing 



^ d i zid i z 2 g{x x ,x 2 ]Zi,z 2 ) a p 1 xG 1 (zi,z 2 ;yi,y 2 )~ 1 \ pa = 5 ap 8 fia 5 A (xi - yi)5 A (x 2 - y 2 ), (3.27) 
one gets 

K(x 1 ,x 2 ;zi,z 2 ) a f3 T s = J d^u^u-Mx (xi , x 2 ; m , u^^xQ^ 1 (ui , u 2 ; z u z 2 )^x T s- (3.28) 

Upon inserting the above expression into Eq. (3.26), we finally arrive at 

K(x 1 ,x 2 ;y 1 ,y 2 ) a p pa = Q(xi,x 2 ;yi,y 2 ) af 3 pa - S(x 1 ,x 2 ;y 1 ,y 2 ) a p pa (3.29) 

where 

S(x ll x 2 ;y ll y 2 ) af3p(T (3 30) 

= J d i u 1 d 4 U2d 4 v 1 d 4 V2Hi(x ll X2;u 1 ,U2) a i3 7 \G~ 1 {ui,U2;v 1 ,V2)- ( \ T s'H 2 (v 1 ,v 2 ;y ll y 2 ) TSpa . [ ' ' 

As we see, the B-S kernel shown above contains two terms. The second term has been explicitly written out. It 
contains a few types of the Green functions as well as the self-energies appearing in the mutually conjugate functions 
7ii{xx,x 2 ; 2/1,2/2) and H 2 (xi,x 2 ;yi,y 2 ) shown in Eqs. (2.35)-(2.40) and (3. 20)-(3.25) respectively. Clearly, to 
give a final expression of the B-S kernel, according to the definition in Eq.(3.2), we need to calculate the function 
Q(xi, X2] 2/1, 2/2) by means of the B-S equations with respect to 2/1 and 2/2 for the Green functions involved in the 
function Hi(xi,x 2 ;yi,y 2 ). We leave the calculations in the next section. 



IV. CLOSED EXPRESSION OF THE B-S KERNEL 



This section is used to calculate the function Q(x\ , x 2 ; yi , y 2 ) so as to give the final expression of the B-S kernel. 
In accordance with the expressions shown in Eqs. (3.2) and (2.35), we can write 

5 

Q(xi,x 2 ;yi,y 2 ) af 3 pa = ^ Qi{xi,x 2 ;yi,y 2 ) a pa (4.1) 

where 

Qi{xi,x 2 ;yi,y 2 ) af 3 pa = [H^\i*d yi +mi - T,)(id V2 +m 2 - iZ c )] a(3p(T (x 1 , x 2 ; 2/1, y 2 )- (4.2) 

This expression indicates that to compute each Qi(xi, x 2 ; 2/1, 2/2)1 we need to derive the B-S equations with respect 
to 2/1 and 2/2 for the Green functions contained in the function H.i(xi,x 2 ;yi,y 2 ). We will separately calculate each 
Qi(xi,x 2 ;yi,y 2 ) below. 

a. Calculation of Qi(xi, X2; yi, y 2 ) 

Based on Eqs. (4.2) and (2.36), we can write 

Qi{xi,x 2 ;yi,y 2 ) af 3 ptT ^ ^ 

= -(r aA % 7 / d 4 u 2 T, c (x2,U2)i3\[G p (^d yi +m 1 - T)(id V2 + m 2 - S c )] 7 a P( j(xi | Xi, u 2 ; 2/1, y 2 ). 

To calculate the above function, it is necessary to start from the following equation of motion which is established 
from Eq. (C6) by setting the source J = 0, 
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G%(xi | xi,u 2 ;yi,y 2 ) 7 \ Ta (i d yi +m 1 ) Tp = -5 lp 5 4 {xi -y 1 )A c p a (x 1 \ u 2 ,y 2 ) Xo - 
-C pa 5 i {y 1 - y 2 )A™(x 1 \ x u u 2 ) 7 a + Gf v {x±; Vl \ x ll u 2 ;y 1 ,y 2 ) lXTa {Y a ^) Tp 

where the Green functions have been represented respectively in Eqs. (2.19), (2.20) and (2.22). Considering the 
definitions given in Eq. (2.26) and in the following 

G${xi\Vi I x\-,u 2 \y u y 2 ) 1 \ pa = Gf v {xi;yj \ x 1 ,u 2 ;y 1 ,y 2 ) 7 x pa 
+A* a (xj | xi,u 2 ) 7 aA* (yj | yi,y 2 ) pa - 

— 1, 2) which immediately follows from Eq. (2.23) and employing Eqs. (3.7) and (4.4), it is not difficult to derive 
[G%(i*d Vl +m l - Y l )] 1 x pa (xi \ xi, u 2 ; yi,y 2 ) = -5~ /p S 4 (x 1 - yi)A c p a (x 1 \ u 2l y 2 ) Xa 



+G p ]l(xi;yi | xi,u 2 ;yi,y 2 ) 7AT(T (r a ^) T(9 - $ d 4 v 1 g c f l (x 1 | xi,u 2 ;vi,y2)-,\T<rV(vi,yi) T p. 

Now, let us operate on the above equation with (i d V2 + m 2 — S c ), giving 

[Gp(i*d Vl +m 1 - Y,){id y2 +m 2 - £ c )] 7A(9(T (xi | xi,u 2 ; yi,y 2 ) 
= -5 lp 5 A {xi -yi)[A p a {id y2 +m 2 - £ c )] A(T (xi | u 2 ,y 2 ) 
+ [Gf L ,(id V2 +m 2 - E c )] 7Ar<T (xi;?/i | xt, u 2 ; y u y 2 ){Y bv ) Tp 
-J d 4 v 1 [g < p \(i d V2 +m 2 - £ c )] 7Arcr (xi | xi, u 2 ; vi, y 2 )T,(vi, yi) rp . 

The first term can be calculated by using the equation derived in Eq. (D5), 

[ A £ a (^s/2 +m 2 - S c )] Act (xi I u 2 ,y 2 ) = A™ b (xi;y 2 | u 2 , y 2 )xs^ hV )sa 
- f d 4 v 2 Af(x 1 | u 2 ,v 2 )\ S Y, c (v 2 ,y 2 ) S a 



(4.6) 



(4.7) 



(4.8) 



where the two Green functions were represented in Eq. (2.19). The third term in Eq. (4.7) can be calculated by 
virtue of the equation given in Eq. (D10) 

G a p {xi | x 1 u 2 ;vi,y 2 )~ f \ T s(i < dy 2 + m 2 ) 5a = -5 Xa 5 4 (u 2 - y 2 )A*(x 1 \ Xi,Vi)~, T ^ 
-C Ta S 4 (v 1 - y 2 )A* a (xi \ x 1 ,u 2 )^ x + Gf u {xi;y 2 \ x 1 ,u 2 ;v 1 ,y 2 ) lX T&(Y V )sa 

in which the Green functions have been respectively defined in Eqs. (2.18), (2.20) and (2.22). In view of the definitions 
in Eqs. (2.26) and (4.5) and employing the equations in Eqs. (4.9) and (3.11), it can be found 

[G p {id V2 +m 2 - Y, c )] lXTa (xi | x\u 2 ;vi,y 2 ) = -8 Xr j5 4 (u 2 - y 2 )A^(xx | x±,v\) 1T 
+Gll(xi;y 2 \ x 1 ,u 2 ]v 1 ,y 2 ) 1 x T s(T ")sa - J d 4 v 2 g a ^{xi \ x\ u 2 ; vi, v 2 ) 7 \ tS ^ c {v 2 , y 2 )sa- 

For calculating the second term in Eq. (4.7), it is needed to use the following equation obtained from Eq. (Dll) by 
setting J = 



[G™(id y2 +m 2 )] 7 xra(xi;yi \ xi, u 2 ; yi, y 2 ) 7ArCT = ~5 Xcr S 4 (u 2 - y 2 )A ab v {xi; y t \ xi,yi) JT 
-C Ta 5 4 (y 1 - y 2 )A™ b (x 1 ;y 1 \ xt, u 2 ) lX + Gf v ° K (xr, yi, y 2 | x 1 , u 2 ; y u y 2 ) 7 \ T 5(^)6,7 



(4.11) 



where all the Green functions were defined respectively in Eqs. (2.18), (2.20) and (2.22). By considering the definitions 
given in Eq. (4.5) and in the following 

Gf v c K (xr,yi,y2 \ x 1 ,u 2 ;y 1 ,y 2 ) lX TS = Gf„ c K (xi;yi, y 2 \ x 1 ,u 2 ;y 1 ,y 2 ) jXT s 
+A* a (xi | xi,u 2 ) 7A A* K c (yi,y 2 | yi,y 2 ) Tr7 

which is implied by Eq. (2.23) and utilizing Eqs. (4.11) and (3.17), an equation of motion for the Green function 
G p t(xi;yi | x\,u 2 ;y\,y 2 ) will be found to be 

[G^i{i*9y 2 +m 2 - £ c )] 7Arcr (xi;?/i | Xi,u 2 ;yi,y 2 ) 7 \ Ta - = -5 Xa 5 4 (u 2 - y 2 )A ab v (x±;y± | x\,y{) 1T 

-C Ta S 4 (y 1 -y 2 )A*^ / b (x 1 ;y 1 \ x lt u 2 ) 7A + G$? K (xi; yi, y 2 \ x 1 ,u 2 ;y 1 ,y 2 ) lXTS {T CK ) Sa (4.13) 

-/ d 4 v 2 Qf v {x\\yi | x\ .u 2 ;y 1 ,v 2 ) 1 x T sT, c (v 2 ,y 2 )sa- 
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When Eqs. (4.8), (4.10) and (4.13) are inserted into Eq. (4.7) and then Eq. (4.7) is inserted into Eq. (4.3), an 
explicit expression of the function Qi(x\, x 2 ;y\, y 2 ) will be derived as given in the following 

Qi(xi,x2\yi,y2) a 0po = Ri(xi,x2;yi,y2) a ppo + M 1 (x i,x 2 ; yi,y2) a P a , . 

+ / d 4 v 1 d 4 v 2 n[ 1 \x ll x 2 ;vi 1 v 2 ) a/}TS 'E(v ll y 1 ) Tp T l c (v 2 ,y 2 )sa- 

where 

Tli{x u x 2 ]yi,y 2 ) a p pa = S 4 (xi - yi){T a ^) ap J d 4 u 2 Y, c (x 2 ,u 2 )p X [A c p a J'(xi;y 2 \ u 2 ,y 2 ) X s(T b ")sa 
-Jd 4 v 2 Al a { Xl \u 2l v 2 ) u iZ c {v 2l y 2 ) Sa ] + iZ c {x 2 ,y 2 ) (ja {T^) ai [Kf u (x l] y 1 \x u yx)- 1 r{T bu )rp ( 415 ) 
-/ d 4 v 1 K p \{x 1 I x 1 ,v 1 ) 1T T,{v 1 ,y 1 ) Tp ] 

S 4 (yi - y 2 )(T a ^) aj J d i u 2 zZ c (x 2 ,u 2 ) fjx A;f(x 1 ;y 1 \ x u «2) 7 a(CT 6 ") <t „ ) 



M 1 {x 1 ,x 2 ;y 1 ,y 2 ) af!pa = -(T afl ) ai J d 4 u 2 Yf{x 2l u 2 ) p\[Gfy c K {xr, yi, y 2 \ x 1 ,u 2 ;y 1 ,y 2 ) 7XrS (T b,/ ) Tp (T CK )sa 

-J d 4 v 2 G ab v ix x -y x I xi u 2] y 1 ,v 2 ) yXrS iZ c (v 2 ,y 2 )sa(T bl ') T p (4.16) 

- f d 4 v 1 gf l ,(x 1 ;y 2 \ x 1 u 2 ;v 1 ,y 2 )~ fXTS Z(v 1 ,y 1 ) Tp (r U )sa] 

and the function Ti^ 1 {x\,x 2 \v\,v 2 ) was defined in Eq. (2.36). 
b. Calculation of Q 2 (x\,x 2 ;yi,y 2 ) 
From Eqs. (4.2) and (2.37), one can write 

Q 2 {xux 2 ;y\,y 2 ) a 0pa ^_ ^_ UVJ) 

= -(r aAI )/3A / d 4 itiS(a;i,tii) Q7 [^(^ yi +m 1 - V)(id y2 +m 2 - Y> c )] lXpa {x 2 | ui, x 2 ; 2/1, 2/2) 

In comparison of the above expression with that denoted in Eq. (4.3), it is clearly seen that the function 
Q 2 (xi,x 2 ;yi,y 2 ) a pp rT can explicitly be written out from the expression of Qi(xi,x 2 ]yi,y 2 ) a p pcr shown in Eqs. 
(4.14)-(4.16) by the following replacements: {T afJ -) aj J d 4 u 2 T< c (x 2 ,u 2 )p X — > (r )p\ J d i u\H{x\, U\) ai and in the 
[G%(id yi + m\ — T,)(id y2 + m 2 - T 1 c )] lXpcr (xi \ X\,u 2 ;yi,y 2 ) whose expression is given in Eqs. (4.7), (4.8), (4.10) 
and (4.13), the gluon field coordinate x\ is changed to x 2 and the quark and antiquark coordinates x\ and u 2 are 
respectively changed to u\ and x 2 . The result is 

Q2(xi,x 2 ;yi,y 2 ) af 3 pa = K 2 (x 1 ,x 2 ;y 1 ,y 2 ) al3pa + M 2 (x 1 ,x 2 ;y 1 ,y 2 ) al3pa 
+ J d 4 vid 4 v 2 H^\xi,x 2 ]Vi,v 2 ) a p T $iZ(vi,yi) Tp Z c {v 2 ,y 2 )8a 



where 



K 2 {xi, x 2 ; yi, y 2 ) a f3 P a = 5 4 {x 2 - y 2 )(T ) /3<7 J d 4 u 1 T.(x 1 , ui) Q7 [A^(x 2 ; j/i \ u ll y 1 )~ fT (T b ") Tp 
- f d 4 v 1 A%(x 2 \u ll v 1 ) 1T ^{v ll y 1 ) Tp ] + i^{x 1 ,y 1 ) ap ^) l 3x[A p \ ab {x 2 \y 2 \ x 2 , y 2 )\s(T V ) S a (41Q] 
-J d 4 v 2 A^ a (x 2 I x 2 ,v 2 )\s'£ c (v 2 ,y 2 )sa} V ' ' 

S 4 (yi -y 2 ){T atl )f }x jd 4 u 1 ^{x 1 ,u 1 ) ai k*- b {x 2 ;y 1 \ u u x i ) lX {CT bv ) ap , 



M 2 (xi,x 2 ;y l7 y 2 ) a0pcr 

= -(r aM )/3A / ' d 4 u 1 H{x 1 ,u 1 ) ai [gf v ° K {x 2 ;y 1 ,y 2 \ u 1 ,x 2 ;y 1 ,y 2 ) jXrS (T bu ) rp (T CK )sa 
- j d 4 v 2 G ab v {x 2 ;y 1 \ u 1 x 2 ]y 1 ,v 2 ) lXT s^ c (v 2 ,y 2 )s r r(T bl/ ) Tp 
-J d 4 viG ab v {x 2 ;y 2 \ mx 2 ; v lt y 2 )jX T S^{vi, yi) Tp (F U )sa] 

( 2 ) 

and the function Tt\ (x\,x 2 ;vi,v 2 ) was defined in Eq. (2.37). 
c. Calculation of Q^{xi, x 2 ; yi, y 2 ) 

In accordance with Eqs. (4.2) and (2.38), the function Qs{xi-,x 2 ;yi,y 2 ) is known to be 

Qz{xi,x2;yi,y2)appa 

= {^) a7 (^ aV )0x[G p t( i ' t 9yi +mi- Z)(id V2 +m 2 - T. c )] yXp a(xi, x 2 \ x u x 2 ; yil y 2 ) 
For calculating this function, it is necessary to start from the equation of motion given by Eq. (C7) 
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G^ b v (xi,X2 \x 1 x 2 ;y 1 ,y 2 ) 7 XTa{id yi +m 1 ) Tp = -5 lp 5 4 (x 1 -y 1 )A c p a J>(x 1 ,x 2 \x 2 ,y 2 )\a (4 22) 

-C pa S 4 (yi -y 2 )A* ll f{x 1 ,x 2 \ x u x 2 ) lX + G ab ^ K {x u x 2 ; y x \ x u x 2 ; y u y 2 ) 7 \ T ^ )r P 

where the Green functions were respectively represented in Eqs. (2.19), (2.20) and (2. 22). Grounded on the definition 
given in Eqs. (2.31) and the following one which can be written out from Eq. (2.23): 



G^k(xi,X2,Ui I xx,x 2 ;yi,y 2 )^x Ta - = G^ b u c K (xi, x 2 , y t \ x 1 ,x 2 \yi,y2) 1 \ra 

+JA* p b { x l, x 2 I X 1 ,X 2 ) 1 xK C (Vi I yi,V2)ra 

here i = 1, 2 and the equations shown in Eqs. (4.22) and (3.7), it is easy to derive 

[S p i(id yi +m 1 - T,)]jx P a{xi,X2 | Xix 2 ;yi,y 2 ) = -5 lp S i (x 1 - y 1 )A c t ^{x ll x 2 | x 2 ,y 2 )xa 
+S P t c K ( x i,x 2 ;yi | x 1 ,x2;yi,y2) 7 \Ta(r CK ) T p - J d 4 v\Q ab v (x\,x 2 \ x 1 ,X2;vi,y2)- r \ T <rZ(v 1 ,y 1 ) 7 

Now, we operate on the above equation with (i d V2 + m 2 — £ c ), giving 

[Gfy(id yi +mi - Y,)(id y2 +m 2 - S c )] 7Ap(T (xi, x 2 \ x 1 x 2 ;y 1 ,y 2 ) 
= -5 7p S 4 (x 1 - y 1 )A^l b (id V2 + m 2 - Y, c )] x<y (xi,x 2 \ x 2 ,y 2 ) 
+ [%l c K { i d y2 +m 2 - T, c )] lX Ta{x u x 2l y 1 \ x u x 2l y u y 2 ){T CK ) Tp 
-J d 4 v 1 [gf l/ (i d V2 +m 2 - T, c )} 7 xra{xi,x 2 | x 1 ,x 2 ;v 1 ,y 2 )T,(v 1 ,y 1 ) Tp . 



(4.23) 



(4.24) 



(4.25) 



The first term is easily computed by virtue of the equation given in Eq. (D6) by adding the term associated to the 
self-energy 

A ™ b (dy2 +™2 - ?, c )\<y{x ll x 2 \ x 2 ,y 2 ) = -5 X(T 5 4 (x 2 - y 2 )i£\ ab u {x 1 - x 2 ) 
+A c p a J'c(x u x 2 ,y 2 \ x 2} y 2 )xs(T CK )sa - J d 4 v 2 A c p a u b (x 1 ,X2 \ x 2 , v 2 ) xs T, c (v 2 - y 2 ) Sa 

where 

iA${ Xl - X2) = (0+ \T[Al{ Xl )Al{x2)]\ 0") (4.27) 

is the gluon propagator and the other Green functions were defined in Eq. (2.19). The third term in Eq. (4.25) will 
be evaluated with the aid of the equation given by Eq. (D12) 



G ab v (x 1 ,x 2 | xix 2 ;v 1 ,y 2 ) 1 xrs(id y2 +m 2 ) Sa = -5 Xa 5 4 (x 2 - y 2 )A ab v (x 1 , x 2 | xi,t>i) 7T 
-CreS 4 ^ - y 2 )A™ b (x 1 ,x 2 I x 1 ,x 2 ) J x + G p [l c K (x 1 ,x 2 ;y 2 \ xi,x 2 ;v 1 ,y 2 ) 1 x T &{T C ' i )sa- 



(4.28) 

Noticing the definitions in Eqs. (2.31) and (4.23) and applying the equations in Eqs. (3.11) and (4.28), wc find 

(4.29) 



d y2 +m 2 - T, c )] 1 xtc{xi,x 2 \ xix 2 ;vi,y 2 ) = -5\ a 5 4 (x 2 - y 2 )A ab u {x u x 2 \ xi, t>i)- 



J-fT 

+Gfv K {xl,X 2 ;y 2 | X l7 X 2 ;V 1} y 2 ) 7 Xrs{T CK )sa - / d 4 V 2 g p \l{x U X 2 I X 1 ,X 2 ;V 1 ,V 2 ) 7 XrS^ C (v 2 ,y 2 )s ! y 

The second term in Eq. (4.25) can be calculated by means of the equation derived in Eq. (D13) 

Ga b u K {xi,x 2 ;yi \ xix 2 ;y ll y 2 )~ t XT8{i*dy 2 +m 2 ) Sa = -8x a S 4 {x 2 - y 2 )A'± c K (xi, x 2 ; yi \ Xi,yi) 7T 
-C Ta S 4 (y 1 - y 2 )A* p a J>c(x 1 ,x 2 ;y 1 \ x u x 2 ) jX + G ab ^ g (x u x 2 ; y x , y 2 \ x 1 ,x 2 ;yi,y2) 7 Xrs(T )sa 



(4.30) 



where all the Green functions can be read off from Eq. (2.18), (2.20) and (2.22). Considering the definitions given in 
Eq. (4.23) and in the following 



Sfj^eixi ,X2;yi, y2 I x u x 2 ;y 1 ,y 2 ) 1 x Ta = G ab ^ g {x 1 ,x 2 ;y 1 ,y 2 \ xi 1 x 2 - 1 y 1 ,y 2 ) lX7 
+^l b {xi,x 2 I x^x^xAZ (y u y 2 \ yi,y 2 )ra 



(4.31) 



which follows from Eq. (2.23) and applying the equations in Eqs. (4.30) and (3.17), it is easy to find 

[G$k(i*dyi +m 2 -Y> c )]~ f x Ta {xi,x 2 ]yi \ xix 2 ;yi,y 2 ) = -5x<?5 4 (x 2 - y 2 )Af u c K (xi, x 2 ; y 1 \ Xi,yi) 1T 

-C Ta 5 4 (y 1 - y 2 )A*^ c (xi,a;2;yi | x 1 ,x 2 ) 1 \ + G^ e {xi,x 2 ;y\,y2 \ xi, x 2 ; yi, y2) 1 \rs{Y )s a ( 4 - 32 ) 

-/ d 4 v 2 gf v c K (xx,x 2 \yx | xi,a;2;yi,W2)7Ar5S c (w2,y2)5<T. 
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Upon substituting Eqs. (4.26), (4.29) and (4.32) in Eq. (4.25) and then inserting Eq. (4.25) into Eq. (4.21), we 
eventually arrive at 

Qz{xi,x 2 - 1 yi,y2)afjpa = K^\x ll x 2 ;yi,y2) a /3pa + K(x 1 ,x 2 ;yi,y2)ap P a + ^3(^1, x 2 ;yi, y 2 ) a ppo- ^ 
+M 3 (x 1 ,x 2 ;y 1 , y 2 ) a p P a + J d i v 1 d i v 2 Hf\x 1 ,x 2 ;vi,v 2 ) al3TS Y>(v 1 ,y l ) Tp Y< c (v 2l y 2 ) Sa 

where 

Ki°\x u x 2 ; yi ,y 2 ) a0pa = 6 4 (xi - yi)<5 4 (x 2 - y 2 )(r^) Q ,(r a V^^i - x 2 ) (4.34) 
which is the t-channel one-gluon exchange kernel, 

^(zi.a^yi.lfiOa/S/ar = {T a ») ai (T av )f 3 xQf v c * e {x 1 ,x 2 -,y 1 ,y 2 \ x 1 ,x 2 ;y 1 ,y 2 ) lXTS {T CK ) Tp {T de )sa, (4.35) 

K 3 (x 1 ,x 2 ]y 1 ,y 2 ) al3prT = -S 4 (x 1 - y 1 )(T a ^) ap (T a ^)p X [A c p a J'^(x 1 ,x 2 ;y 2 \ x 2 ,y 2 ) xs {T cli )sa 
-J d 4 v 2 A™ b (xi,x 2 I x 2 , v 2 )\sY. c (v 2 , y 2 )sa] 

-5\x 2 ~y 2 ){t v ) Pa ({T^) ai [kf u c K (x ll x 2 -y 1 \ x u yi) 1T {T™) Tp (4.36) 
-/ d 4 v 1 A p [ b / (x l7 x 2 I a;i,Ui) 7T S(t;i,j/i) r p] 

+5\yi ~y 2 ){T^) ai {T au )p X Alfc{x u x 2 -, yi \ x u x 2 ) 7X (CT CK ) ap , 



M3(x 1 ,x 2 ;yi 1 y2)a0pa 

= -(r o ")a 7 (r° , ') / 3A[/ d A v 2 g a ^(x 1 ,x 2 -y 1 | ii,^;?/!,^^^,^)^^)^ (4.37) 
+ f d 4 v 1 g p lb l/ c K (x 1 ,x 2 ;y 2 | xi, x 2 ; vi, y 2 ) 1 \r6^(vi, yi) Tp (r CK )sa} 

and the function H^\xi, x 2 ; v\, v 2 ) was defined in Eq. (2.38). 
d. Calculation of Q&{x\,x 2 \y-i,y 2 ) 

According to Eqs. (4.2) and (2.39), the function Qi{x\,x 2 ;y\,y 2 ) is ought to be 

Qiixi.x^y^y-^jappa _^ ^_ ^_ , 43g , 

= <5 4 (xi - x 2 ){T ail C){3 a [A* p a {id yi +m 1 - iZ)(id V2 + m 2 - T,°)] pr7 (x 2 | 2/1,2/2) 

Clearly, to evaluate the above function, we need at first the equation of motion which is obtained from Eq. (Cll) by 
setting the source J to be zero and adding the term associated with the self-energy 



[A* a («3yi - Z)] pa (x 2 1 2/1,2/2) = A*7(z2;yi 2/1, y2)ra{T hv ) Tp 

-J d A Vik* p (x 2 I Ul,2/2)raS(wi, 2/1)7 



b, 



(4.39) 



l-rp 



where the Green functions was written in Eq. (2.21). 

Now, let us operate on Eq. (4.39) with the operator (i d V2 + m 2 — S c ), 

[K* p a {id Vl + mi - Y){id V2 +m 2 - T. c )] pa (x 2 | 2/1,2/2) 

= [A;; b ^ 2 + m 2 - £ C )] TCT (* 2 ; 2/1 I 2/1, V2)(T*»')t P (4-40) 
- / d 4 v 1 [A p a (i d V2 +m 2 - S c )] T(T (x 2 | v 1 , ?/ 2 )£(ui, yi) Tp . 

Based on the equations shown in Eqs. (D18) and (D19), one can write respectively the following equations: 



[^{idy 1 + m 2 - £ c )] T(T (x 2 I v 1 ,y 2 ) = A*" b (x 2 ;y 2 | «i, 2/2 )r«(r "),$,, 
-/ d 4 u 2 A* a (x 2 I ui,w 2 )r^S c (w 2 ,2/ 2 )fc 



and 



r -r-*a6 , 



[ A ^7(*^y 2 + TO 2 - S c )] TCT (x 2 ,2/i I 2/1,2/2) = C TCT (5 4 (2/i - y 2 )iA^,(x 2 - 2/1) 
+^"^(^2,2/1,2/2 I 2/i,2/2)T5(r CK ) 5( T - fd 4 v 2 A* l 2 (x 2 ,yi I 2/i^2)r«5S c (w 2 ,y 2 ) 5o 



(4.41) 



(4.42) 



where iA p 'l(x 2 — 2/1) is the gluon propagator as defined in Eq. (4.27) and the other Green functions can be read from 
Eq. (2.21). 
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Substituting at first Eqs. (4.41) and (4.42) into Eq. (4.40) and then Eq. (4.40) into Eq. (4.38), the function 
Qi{x\,x 2 \ 2/1,2/2) will be explicitly written in the form 

Qi{xi,X2\yi,y2)cti3pa = Ki°\xi,x 2 ;yi,y2) a f3po- + K i (xi,x 2 ;yi 1 y2) a i3pa ^ ^ 

+ J d i v 1 d i v2Ti ( i\xi 1 x 2 ;vi 1 v2) a f3T8'E(vi 1 yi)Tp?> c (v2 1 y2) 

where 

Kf\x u x2\yi 1 y2) a ppa^5\x 1 ~x2)5\y 1 -y2)(T a ^C) Pa {C- 1 T bv ) ap i^f l/ {x2-y 1 ) {AAA) 
which is the s-channel one-gluon exchange kernel, 

T^i(xi,x 2 ; 2/1, y 2 ) a f3pa = S i (x 1 - x 2 ) (T aft C)p a { [A*" fc K c (x 2 ; 2/1 , 2/2 I Vi, ^riP^ta 

-JdS 2 A*; h (x 2 ;yi I 2/i,«2)r^ c («2,2/2W](r bl rp (4.45) 

-JdVA*" (x 2 ;2/ 2 I v 1 ,y 2 ) T s^(v 1 ,y 1 ) Tp (T ' / ) Sa } 

and the function Ti.^ (x\, x 2 ; v\, V2) was defined in Eq. (2.39).. It is noted here that the s-channcl one-gluon exchange 
kernel describes the qq annihilation and creation process which takes place between the quark (antiquark) in the initial 
state and the antiquark (quark) in the final state as indicated by the gluon propagator in Eq. (4.44), not between the 
quark and the antiquark both of which belong to the initial state or the final state B-S amplitude for a bound state. 

e. Calculation of Qs(xi, X2; J/i, 2/2) 

From Eqs. (4.2) and (2.40), we have 



Q5(xi,X2;yi,y2)af3p<T 

= J d 4 u 1 d 4 u 2 T,(x 1 ,u 1 ) ai Y 1 c (x 2 ,u 2 )i3\[Q(i d Vl + m x - d V2 + m 2 - T, c )] 1 x pa {u ll u 2 ; 2/1, 2/2)- 
By applying the equation given in Eq. (3.19), the above function can explicitly be represented as 



(4.46) 



Q5(xi,x 2 ;yi,y 2 ) a/ 3p a = U (xi,x 2 ;yi,y2) a f3pa , 4 ^ 
+ J d 4 uid 4 u 2 Y,(xi,ui) a ^ c (x 2 ,u 2 )0\TC 2 (ui,u 2 ;yi,y 2 ) 7 \ pa 

where 

T^o(xi,x 2 ;y l7 y 2 ) a p pa = Y:(x 1 ,y 1 ) ap j: c (x 2 ,y 2 )p a (4.48) 

and H 2 (u\,u 2 ; 2/1, y 2 )-y\ pa was shown in Eqs. (3.20)-(3.25). 

In summary, combining the expressions given in Eqs. (4.14), (4.18), (4.33), (4.43) and (4.47), an explicit expression 
of the function Q(xi,x 2 ; 2/1,2/2) written in Eq. (4.1) will be eventually found as shown in the following 



Q{x-i,X2-,yi,y2) a ppa = K^\x 1 ,X 2 ;yi,y2)al3pa +K { S '?\xi,X 2 ;yi,y2)al3p<7 

+K(X1 , X2\ 2/1 , 2/2)a/3p<T + H(X\ , X2 \ 2/1 , 2/2)a/3 pC r + M (xi , X 2 ] 2/1 , 2/2)a/3p CT 

+ / d 4 u 1 d i u 2 T,(x 1 ,u 1 ) a7 J: c (x2,u2)pxH2(u 1 ,u2;yi,y2) 1 xpa 

+ J d 4 V 1 d 4 V 2 H 1 (x 1 ,X2;V 1 ,V2) al 3rS^(vi,yi)rp^' C (v2,y2)sa 



(4.49) 



where we have defined 

4 

R-(xi,x 2 ;yi,y2) a f3pa = ^'R- l {xi,x2;yi,y2) a t3pa, (4.50) 



i=0 



M(x 1 ,X2;yi,y 2 ) a ppa = ^2M i (x 1 ,x 2 ;yi,y2) a fipa (4.51) 

and 

H j (x 1 ,x 2 ;yi,y2) a ppa = {x u x 2 - 1 yi,y 2 ) afjpa (4.52) 



i=i 



4 

AO, 
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here j = 1,2 and all the functions on the RHS of Eqs.(4.50)-(4.52) have already been explicitly given before. 

Now let us turn to the function <S(a;i, X2\ J/i, V2)ai3po expressed in Eq. (3.30). According to the definitions denoted 
in Eqs. (2.35), (3.20) and (4.52), we can write 

Hj{x 1 ,X2;yi,y2) a pp<7 = Hj(xi, x 2 ; yi, y 2 ) a fi P a + Tif\x u x 2l j/i, y 2 ) a ppo (4.53) 
.Thus, Eq. (3.30) can be rewritten in the form 

S(xi,X 2 ;yi,y2)afjpcr 

= j d i uid i u 2 d i v 1 d i V2H 1 (xi,x 2 ;ui,U2) a i3- f \G~ 1 (ui,u 2 ;vi,v 2 )~ f \ TS H2(vi,V2;yi,y2)TSpa 
+ j d i u 1 d i u 2 d i v 1 d i v 2 n[ 5 \x 1 ,x 2 ;u 1 ,u 2 ) a ^xG~ 1 (ui,u 2 ;v 1 ,v 2 )^x T5 n 2 (vi,v 2 ;yi,y 2 ) T s P a 
+ J d\ 1 d i u 2 d i v 1 d i v 2 H 1 (x 1 ,x 2 ;u 1 ,u 2 ) al 3 1 xg- 1 (ui,u 2 ;v 1 ,v 2 )^x TS n 2 \v 1 ,v 2 ;y 1 ,y 2 ) TSp „. 

When the functions H^\xi, x 2 ; tii, U2) a p-y\ and Ti-^ (vi,v 2 ;yi,y 2 ) T s pa in the second and third terms are respectively 
replaced by the expressions given in Eqs. (2.40) and (3.25) and utilizing the relation in Eq. (3.27), it is easy to see 
that the last two terms in Eq. (4.54) exactly equal to the last two terms in Eq. (4.49) respectively. Therefore, when 
we substitute Eqs. (4.49) and (4.54) into Eq. (3.29), the above-mentioned terms are cancelled with each other. Thus, 
the B-S kernel will be finally expressed by 



K(xi,x 2 ;yi,y2) af>i>a = Q{xi,x 2 ;yi,y2)ai3p^ - S(xi,X2;yi,y2) a 0p<j (4-55) 



where 



Q(xi,x 2 ;yi,y 2 ) afjpa = K^\xi, x 2 ;yi,y 2 ) a f3 P a + Ki" ) (x 1 ,x 2 ;y 1 ,y 2 ) afjpiJ ^ ^ 

+K{x 1 ,x 2 :yi,y 2 ) + 1Z{x l ,x 2 ]y ll y 2 ) afjpa + M(xi, x 2 ; yi, y2)ap P a 

S(xi,x 2 ;yi 1 y 2 ) a ppa (4 57) 

= j d i uid i u 2 d i vid i v 2 H 1 (x 1 ,x 2 \u 1 ,u 2 ) al 3 1 \Q^ 1 {u 1 ,u 2 \v 1 ,v 2 ) 1 x T sH2{v 1 ,V2]yi,y2)rS pa - 

In Eqs. (4.56) and (4.57), all the functions are represented through only a few types of Green functions. Therefore, the 
kernel derived is not only easily calculated by the perturbation method, but also provides a proper basis for studying 
the QCD nonperturbative effect. 



and 



V. SUMMARY AND DISCUSSIONS 



In the preceding sections, the B-S kernel has been derived by means of the B-S equations satisfied by the qq 
Green functions and given a closed and explicit expression shown in Eqs. (4.55)-(4.57). This expression is manifestly 
Lorentz-covariant even in the space-like region of Minkowski space where the bound states exist and it suits to not 
only the case that the quark and antiquark have different flavors, but also the case where they have the same flavor. 
In the latter case, the masses of the quark and antiquark are the same, mi = m 2 = m. Meanwhile, the B-S kernel 
contains the part describing the qq annihilation and creation process such as the s-channel one-gluon exchange kernel 
Ki°\xi,x 2 ;yi,y 2 ) shown in Eq (4.44) and the terms involving the Green functions A*°;;^;_" (a;j , yj... \ xi,x 2 ) and 

^u...v... ( x i---iVj,— I 2/1)2/2) which give the higher order perturbative corrections to the vertices of the s-channel one- 
gluon exchange kernel. It is well known that the s-channel one-gluon exchange kernel gives no contribution to the 
qq bound state because the matrix element of the color operator in the kernel vanishes in the qq color singlet. But, 
this does not imply that the qq annihilation part of the kernel is meaningless. For example, the s-channel one-gluon 
exchange kernel derived in this paper would give a considerable contribution to the interaction of a multi-quark- 
antiquark system as shown in our previous works [37]. For the former case that the quark and antiquark are of 
different flavors, certainly, the qq annihilation part of the kernel is absent. In this case, the B-S kernel is of a simpler 
expression which can be written out from the corresponding expression given in the preceding section by deleting out 
the qq annihilation terms and replacing all the Green functions by the conventional ones. For convenience of later 
discussions, we summarize the expression of such a kernel in the following. 

K(xi,X2;yi,y2) a p P a = Q(x\,x 2 \y\,y 2 ) a(ipa - S{x 1 ,x 2 - 1 yi,y 2 ) a p pcr (5.1) 

where Q(x\, x 2 ; yi, 2/2) and S(x\,xi\yi,y?>) are separately described below. 
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Q{xi,x 2 ;yi,y2)a0 P <T = K\ \x\,x 2 \y\,y 2 ) aS j pa + R(xi, x 2l yi, y 2 ) a p P a ^ 2 ) 

+K(xi,x 2 ;yi, y 2 ) a fi P a + M(x 1 ,x 2 ;y 1 , y 2 ) a (3 P a 

where Kf\xi,x 2 - ) y ll y 2 ) al3pa was given in Eq. (4.34), 

K(x 1 ,x 2 ;y 1 ,y 2 ) a p pa = (T a ^) al (T a ^) 0X G a p b u c ^ (x 1 ,x 2 ;y 1 ,y 2 \ x 1 ,x 2 ;y 1 ,y 2 ) lXT s(T CK ) Tp (T d0 )sa (5.3) 

which formally is the same as given in Eq. (4.35), but the Green function in it becomes the conventional one as 
defined in Eq. (2.22), 

3 

R(x 1 ,x 2 ;y 1 ,y 2 ) a i3pa = ^Ri{xi,x 2 ;y 1 ,y 2 ) aPpa (5.4) 

here the fact that 1Z±(x\, x 2 \ yi,y 2 ) defined in Eq.(4.45) vanishes in this case has been considered, Ro(x\, x 2 ; yi,y 2 ) a {3 P o 
is defined as the same as written in Eq. (4.48), while the other Ri(x\, x 2 ; yi, y 2 ) a i3pa are now defined by 

Ri{xi,x 2 ;y 1 ,y 2 ) al 3 pa = 8 4 {x l - yi){T afl ) ap J d 4 u 2 Y; c (x 2 ,u 2 ) l 3 X {Ap'(x 1 ;y 2 \ u 2 ,y 2 ) X s(T u ) Sa 

- J d i v 2 Kl a {x l | u 2 ,v 2 ) X5 Y, c {v 2 ,y 2 ) 5a ] +^%x 2 ,y 2 )^(T^) a , ( [Af u (x 1 ;y 1 \ x uyi )^{T b ») Tp (5.5) 

- Jd^A^xx | x 1 ,v 1 ) 1T T,(v 1 ,y 1 ) Tp ], 



R2{xi,x 2 ;y 1 ,y 2 ) aj3p<7 = S 4 (x 2 - y 2 ){T ail ) f3<y J d 4 u 1 Y,(x 1 ,u 1 ) aj [Af u (x 2 ;y 1 \ u 1 ,y 1 )~ tT (T bv ) Tp 

-JdV;(i 2 |«i,«i) T rS(wi,yi)rp]+S(a:i,yi) a ^(r '')^ A [A^(a:2;ite \ x 2 ,y 2 ) xs {Y bv ) Sa (5.6) 
-/ d 4 v 2 A c p a (x 2 | x 2 ,v 2 ) xs T, c (v 2 ,y 2 ) Sr7 ] 



and 



Rsixi.x^y^y^a^p^ = -5 4 (x 1 - y 1 )(T a ^) ap (T av ) l3X [A c ^(x 1 , x 2 ; y 2 | x 2 , y 2 ) xs (T CK ) S<7 
- J d 4 v 2 A^f(x 1 ,x 2 | x 2 , v 2 ) x& Y, c (v 2 , y 2 )sa] 



-6 4 (x 2 -y 2 )(T )^{{T^) ai [kf^{x u x 2 - yi \ x u y^) 1T {T™) 
- J d 4 v 1 A^ b l/ (x 1 ,x 2 | x 1 ,v l ) 1T Y,(v u y l ) Tp ] 



(5.7) 



Tp 



which are respectively written from Eqs. (4.15), (4.19) and (4.36) by dropping the terms related to the qq annihilation, 
and 

3 

M{x 1 ,x 2 ;y 1 ,y 2 ) aj3prT = Mi (x\,x 2 ',yi, y 2 ) a p prT (5.8) 



here 



M 1 (x ll x 2 ;y 1 ,y 2 ) a p pa _ 

= -(r^a-y / d 4 u 2 ^x 2 ,u 2 )p X [G«i c K (x 1 -,y 1 ,y 2 I x 1 ,u 2 ;y 1 ,y 2 ) 7XT s(T^) Tp (T CK ) Sa 

- f d 4 v 2 G p [ b u (x 1 ;y 1 \ x 1 u 2 ;y 1 ,v 2 ) 7XT s?, c (v 2 ,y 2 )s < j(r bl/ ) T p ^ ' > 

-J d 4 v x G ab {x 1 ;y 2 \ x 1 u 2 ;vi,y 2 ) 1 \ T s'E(vi,yi) Tp (r V )sa], 



and 



M 2 (x ll x 2 ;y 1 ,y 2 ) a fi prT 



; 2/1, 2/2 | u 1 ,x 2 ;y 1 ,y 2 ) lXT s(T b ") Tp (r CK )sa . 
j d 4 v 2 G ab v {x 2 ]y 1 \u 1 x 2 ]yi,v 2 ) lXTS Y 1 c (v 2 ,y 2 ) Sa (T bv ) Tp ^' ' 



{IIS \ 

J d 4 v 1 Gf u (x 2 ;y 2 | Ui,a; 2 ; v lt y 2 ) lXT s^(vi, yi) rp (T v ) Sa ] 



M 3 (x 1 ,x 2 ;yi 1 y 2 ) a f ipcr 

= -(T^) a7 (T av ) !3X [Jd 4 v 2 G^ c K (x 1 ,x 2 ;y 1 \ x u x 2 - Vl , ^ At5 S c (« 2 , y 2 )^(r CK ) r , (5.11) 
+ J d 4 v 1 G p [ b / c K (x 1 ,x 2 ;y 2 \ x 1 ,x 2 ;v 1 ,y 2 ) 1 x TS i:(v 1 ,y 1 ) Tp (T CK ) Sa ] 
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which are written respectively from Eqs. (4.16), (4.20) and (4.37) with the Green functions being replaced by the 
conventional ones. 

The function S(xi,x 2 ;yi,y 2 ) is defined by 

S(xi,x 2 ;yi,y2)ai3pcT (5 12) 

= J d i U 1 d i U 2 d A V 1 d' i V2H 1 (xi 1 X 2 ] Ul,U 2 ) a/ 3 7 xG~ 1 (ul,U2;Vl,V 2 ) 7 XT5H 2 (vl,V 2 ; yi, y2)rSpa 



where 



H j (x 1 ,x 2 ;ui,u 2 } a 7 x = ^2Hy ) (xi,x 2 ;yi,y2) a p P a (5.13) 



i=i 



here j = 1,2, 



and 



H[ 1 \x 1} x2;yi,y2)ai3 P a = -(r aM ) Q7 J d i u 2 Y, c (x 2 ,u 2 ) xG^(xi I x 1 ,u 2 ;yi,y2) 7 \ P <T, (5.14) 

H[ 2 \x 1 ,x 2 ;yi,y2)ap P <7 = -(r^A J d 4 u 1 Y,(x 1 ,zi) ai G b u (x 2 \ ui,x 2 ;yi,y 2 ) 1 x P <7, (5.15) 

H i 1 3 \xi 1 x2;yi 1 y2) a 0pa = (r aAI ) Q7 (r fc ^ ')p\Gf u {xi, x 2 \ xi,x 2 ;yi,y 2 ) 7 x P a, (5.16) 

H^\x 1 ,x 2 - 1 y 1 ,y 2 ) a f 3prT = - J d i v 2 G a ^(y 1 \ x 1 ,x2;yi,v 2 ) a firS^ c (v2,y2)8a(r af ")rp, (5.17) 

H 2 2) (xi,x 2 ;yi,y2)afjp<7 = - J d A viG b v {y 2 \ x 1: x 2 ; v 1: y 2 ) a T 5^(vi, yi) T p(F bl ')5<? (5.18) 

(x 1 ,x 2 ;y 1 ,y 2 ) a0pa - = Gf v {y u y 2 \ x u x 2 ; yi, y2) a f3Ts(T af ") T p(r b ' / ) Sa (5.19) 



which are written from Eqs. (2.36)-(2.38) and (3.21)-(3.23) with replacing all the Green functions by the conventional 
ones. In the above, Ti.^ 4 \x\, x 2 ; yi, y 2 ) = has been considered. 

Now, some discussions on the kernel derived in this paper are in order. For simplicity, we concentrate our attention 
on the expressions shown in Eqs. (5.1)-(5.19). First, we note that the terms in the function M(xi,x 2 ;yi,y 2 ) defined 
in Eqs. (5.8)-(5.11) are all related to the quark and antiquark self-energies S(xi, y\) and Yf(x2 1 y2)- These terms are 
responsible for eliminating the corresponding terms contained in the function K(xi,x 2 ;yi,y 2 ) defined in Eq. (5.3). 
This point would be clearly seen from the one and two-particle irreducible decompositions or the perturbative expan- 
sions of the Green functions occurring in the functions K{x\, x 2 ; yi, y 2 ) and M(xi, x 2 ; y\ 1 y 2 ) ( detailed discussions of 
this point will be presented later) . By the irreducible decomposition of the Green function involved in the function 
K(xi,x 2 ;yi,y 2 ), it would be found that the irreducible diagrams of the function K(xi,x 2 ;yi,y 2 ) can be divided into 
two classes: the first class consists of the diagrams with the self-energies of the so-called external quark and anti- 
quark (which mean the quark and antiquark occurring in the B-S amplitude); the second class contains the diagrams 
without such self-energies. The first class of diagrams are precisely cancelled out by all the diagrams of the functions 
M(xi,x 2 ;yi,y 2 ). Similarly, the self-energy-related terms in the function R(xi,x 2 ;yi,y 2 ) defined in Eqs. (5.4)-(5.7) 
and (4.48) are responsible for cancelling the corresponding terms contained in the two self-energy-irrelevant terms 
in the function R^{x\, x 2 \ yi, y 2 ) shown in Eq. (5.7). From this cancellation, we see that the B-S kernel is, actually, 
irrelevant to the self-energy corrections of the external fermion lines. It only contains the terms consisting of vertices 
and internal propagators. In particular, the function R(x\, x 2 ; yi, y 2 ) gives all the higher order perturbative correc- 
tions to the bare vertices included in the kernel Kj.°\xi,x 2 ;yi 1 y 2 ) so as to make the vertices to be exact. Therefore, 
the first two terms on the RHS of Eq. (5.2) actually give the exact one-gluon exchange kernel. The self-energies of 
the external fermions are included in the full propagators Sf(xi — z\) and Sp{x 2 — z 2 ) in Eq. (2.45). With these full 
propagators, the B-S equation has an advantage that it can conveniently be renormalizcd. This advantage interprets 
why in the derivation of the B-S equation, we operate on the Green functions with the operators (id Xl — mi + S) 
and (id X2 — m 2 + E c ) other than with the operators (id Xl — mi) and (id X2 — m 2 ). Certainly, the B-S equation can be 
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set up by using the operators (id Xl — mi) and (id X2 — m 2 ) to operate on the Green function G{x\ 1 x 2 ; yi, 2/2)- In this 
case, the fermion propagators in Eq. (2.45) become free ones and, meanwhile, the B-S kernel, in the case that the 
quark and antiquark are of different flavors, has an expression which may be written down from Eqs. (5.1)-(5.19) by 
deleting out the self-energy-related terms, as shown in the following 

K(x 1 ,x 2 ;y ll y 2 ) = KP(x 1 ,x 2 ;yi,y_2) + R{xi,x 2 ;y 1 ,y 2 ) / 5 20 ^ 

+G(x 1 ,x 2 ;y 1 ,y 2 \ xi,x 2 ;y±,y 2 ) - S{x±,x 2 ;y±,y 2 ) 

where Kf\xi,x 2 ;yi,y 2 ) is still represented in Eq. (4.34), 

R(xx,x 2 ;yi,y 2 ) = S{x 2 - y 2 )A(xi,x 2 ,yi \ x±,y{) - S(xx - yi)A c (x 1 ,x 2 ;y 2 | x 2 ,y 2 ) (5-21) 

here 

A(x 1 ,x 2 ;y 1 \ x uyi ) = (0+ |T{J(a;i)A(a:2)J(»i)}| 0~> (5.22) 



A c (x 1 ,x 2 ;y 2 \ x 2 ,y 2 ) = (0+ T{A( Xl )3 c (x 2 )J c (y 2 )} 



(5.23) 



in which 



A(x) = T a ^A^(x), A(x) -- 
J c (x) = A~(x)ip c (x),~J(x) 



_ Al(x), 3(x) = A(x)1>(x), 
ip(x)A(x), J (x) — tp (x)A(x) 



(5.24) 



and 



G(xi,x 2 ;yi,y 2 \ xi,x 2 ;yi,y 2 ) 



= 0+ 



T{J(x 1 )J c (x 2 )J( 2 / 1 )J (y 2 )} 



(5.25) 



S(xi,x 2 ;yi,y 2 ) 

= / d i u 1 d 4 u 2 d i v 1 d i v 2 G(x 1 ,x 2 \ x 1 ,x 2 ;u 1 ,u 2 )G- 1 (u ll u 2 ;v 1 ,v 2 )G(y 1 ,y 2 \ v 1} v 2 ; y u y 2 ) 



(5.26) 



here 



G(xi,x 2 I xi,x 2 ;ui,u 2 ) = (0+ T{3(xi)3 c (x 2 )^(ui)i' (u 2 )} 



(5.27) 



G( Vl ,y 2 I v u v 2]yi ,y 2 ) = (0+ ^{^(v^rMJiy^iy^ 0") . (5.28) 

The kernel shown in Eqs. (5.20)-(5.28) formally is simple though, its content actually is more complicated than 
the kernel written in Eqs. (5.1)-(5.19) because it contains the diagrams related to the self-energy corrections of the 
external fcrmions. 

It is well known that the B-S equation is invariant with respect to renormalization. In other words, equation (2.45) 
keeps the same form before and after renormalization. Therefore, the renormalized B-S kernel, for example, in the 
case that the flavors of quark and antiquark are different, is still represented by Eqs. (5.1)-(5.19) provided that all the 
quantities such as the quark and antiquark masses, the coupling constant, the self-energies and the Green functions 
are replaced by the renormalized ones. In such a kernel, the first two terms in Eq. (5.2) give the exact renormalized 
one-gluon exchange kernel. Especially, at the renormalization point, as demonstrated in the renormalization theory [5, 
38], the one gluon exchange kernel is reduced to the form as given in the lowest order approximation of perturbation 
and the function M{x\, x 2 ; yi, y 2 ) vanishes because at the renormalization, according to the renormalization boundary 
conditions, the propagators become the free ones due to vanishing of the self-energies and the vertices become the bare 
ones. Only in this case, the B-S kernel has a simple form as written in Eq. (5.20) with the function R(xi, x 2 ; yi, y 2 ) 
defined in Eqs. (5.21)-(5.23) being absent. 

Next, we would like to mention the role played by the last term in Eq. (5.20). As pointed out in Sect. 2, the 
Green function G{x ll x 2 \ x\,x 2 ;u\,u 2 ) with two-gluon fields at the positions x\ and x 2 is B-S reducible and can be 
represented as 
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G(xi,X2 I xi,X2]Ui,U2) = / d 4 z 1 d 4: z2K(xi,x2;z 1 ,z2)G(z 1 ,z2;ui,u2) 



(5.29) 



where K(x\,X2\z\,Z2) is just the B-S kernel which is generated from the Green function G{x\,X2 \ ari,^; iti, 112)- 
Similarly, the Green function G(yi,y2 \ v\, V2\ yi, 2/2) with two gluon fields at y\ and 2/2 is also B-S reducible and can 
be expressed in the form 

G{yi,y 2 I v 1 ,v 2 ;yi,y2) = J d i z[d 4 z 2 G(v 1 ,v 2 ; z[, z 2 )K(z[, z 2 ;y 1 ,y 2 ). (5.30) 

Substituting Eqs. (5.29) and (5.30) into Eq. (5.26) and using the identity denoted in Eq. (3.27), it is found that 

S(xi,x 2 ;yi,y2) = J d 4 u 1 d 4 u2d 4 vid 4 v2K(x 1 ,x2;u 1 ,u2)G(u 1 ,u2;v 1 ,v2)K(v 1 ,v2]yi 7 y2)- (5.31) 

This expression shows the typical structure of the two-particle reducible part of the B-S kernel. On the other hand, the 
Green function G{x\, X2] 2/1, 2/2 | x\, X2; yi, 2/2) for which at every position a gluon field is nested, in general, can be split 
into a B-S irreducible part Gm(xi, X2] 2/1,2/2 | x\, X2; J/i, 2/2) and a B-S reducible part Gre(xi, X2] yi, 2/2 | x\, X2', J/i, 2/2)- 
The B-S reducible part just equals to the function S(xi, X2\ yi, 2/2)- Therefore, both of the functions S(xi, X2; 2/1, 2/2) 
and Gre{x\, X2\ 2/i, 2/2 | £1, £2; j/i, 2/2) are cancelled with each other in Eq. (5.20). Thus, we have 

K(x 1 ,x 2 ;yi,y2) = K t (x 1 ,x 2 ;yi,y2) + G IR (x 1 ,x 2 ,yi,y2 \ x 1 ,x 2 ;yi,y 2 ) (5.32) 

in which the B-S irreducible part of the Green function represents two and more than two-gluon exchange interactions 
in the sense of perturbation theory. Similar cancellations take place in Eqs. (4.55) and (5.1) where the function 
S(xi, X2\ 2/i ; 2/2) and S(xi, X2; 2/1,2/2) just play the role of eliminating the B-S reducible parts contained in the remaining 
terms in the kernels. 

At last, we would like to stress that unlike the Dyson-Schwinger (D-S) equation [39,40] which contains an infinite 
set of equations, the B-S equation is of a closed form. In particular, the interaction kernel in the equation represents 
all the interactions taking place in the bound states. As shown in the preceding section and this section, the closed 
expression of the B-S kernel contains only a few types of Green functions. These Green functions can easily be 
calculated by the perturbation method. In this kind of calculations, we only need the perturbative expansions of the 
Green functions without concerning the calculation of other more-point Green functions as it is necessary to be done 
for the D-S equation. Especially, the Green functions are possible to be evaluated nonperturbatively as suggested 
by the lattice gauge theory. Therefore, the expression of the kernel given in this paper provides a new formalism for 
exploring the QCD nonperturbative effect existing in the bound state and the quark confinement. For this purpose, 
it is appropriate to start from the kernel represented in this section, particularly, as the first stage, from the kernel 
in Eq. (5.32) in which the last term is given by the difference between the G(#i, £2; Z/i, 2/2 I x±, x 2 ] 2/1, 2/2) m Eq. 
(5.25) and the S(x\, X2; 2/1,2/2) in Eq. (5.26). In the ordinary quark potential model, the second term in Eq. (5.32) is 
usually simulated by a linear potential kernel or some other modified confinement one. Obviously, this simulation is 
oversimplified. To search for a sophisticated confining potential including not only its spatial form, but also its spin 
and color structures, it is necessary to call for a nonperturbative calculation of the the B-S kernel given in this paper. 

The spin and color structures are determined not only by the matrices T afl and T in Eq. (5.25), but also by the Green 
functions in the kernel. For example, the color operator for the quark- antiquark system can generally be decomposed 
into a combination of the two particle color operators 1(g) 1, T a <£> 1, 1 (g> T and T a ® T . The coefficients of the 
combination can be determined in a certain perturbative calculation, but can not be concretely given nonperturbativly 
before a nonperturbative calculation of the Green functions in the kernel is implemented. Since the B-S amplitude in 
Eq. (2.44) is color singlet, containing a color singlet wave function in it, to solve the B-S equation for bound states, 
one needs to evaluate the expectation value of the color operator between the color singlet. Clearly, once the B-S 
kernel could be calculated nonperturbatively, the B-S equation can be solved to give not only the exact masses of qq 
bound states, but also the B-S amplitudes which may be used to investigate the meson decay processes. This just is 
the advantage of the formalism of B-S equation. At last, it is mentioned that the procedure of deriving the qq bound 
state B-S kernel can straightforwardly be applied to derive the B-S kernels for other multi-quark and/or multi-gluon 
bound systems as well as other fermion and/or boson systems either at zero-temperature or at finite temperature. 
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VII. APPENDIX: EQUATIONS OF MOTION SATISFIED BY GREEN FUNCTIONS 



In the previous sections, the B-S equations used to derive the B-S kernel were derived based on the equations of 
motion satisfied by some kinds of Green functions. In this appendix, we show how the equations of motion which 
describe variations of the Green functions with the coordinates x\, X2, yi and yi are derived from the QCD generating 
functional. 



A. Equations of motion with respect to xi 



First, wc write down the QCD generating functional [5] 

Z[J,1f,V,Zt] = ]~ j V(A,^^C,C)e u (Al) 

where 

1 = J d 4 x[C + J a »A; + rft> + + lC + C$ (A2) 
in which C is the effective Lagrangian of QCD 



C = 4>(id-m + gA)ip - - F a » v F* v - -L(0MJ1) 2 + C a d^{D ab C b ) (A3) 



here A =j^T a A^ with A^ being the vector potentials of gluon fields, 



F% = d»K - d v A« + gf ahc A\Al (A4) 

are the strength tensors of the field, 

Df = 5 ab d„ + gf abc AlC c (A5) 

are the covariant derivatives, C a ,C b are the ghost fields, and J aAt ,ry, 77, £ and £ denote the external sources for the 
gluon, quark and ghost fields respectively. By the charge conjugation transformations shown in Eq. (2.2) for the 
quark fields and in the following for the external sources 

rf = CrJ T ,rf = -rfC' 1 , (A6) 

it is easy to prove the following relation 

ip(id—m + gA)ip + rjip + tpr) = ip c (id-m + gA)ip c + rfip c + ^ c r\ c (A7) 

where A = 7^T°A" with T° = — A a */2 being the antiquark color matrix. 

Now let us proceed to derive the equations of motion obeyed by the Green functions. Taking the functional derivative 
of the generating functional in Eq. (Al) with respect to the field function i\) a {x\) and considering 

5Z = (A8) 



and 

SI 



we have 



= Va(xi) + [{id Xl - mi) Q7 + 3A(xi) Q7 ]V'7(xi), (A9) 



{Va ( Xl ) + [{ i8 Xl - m lU + (^^^O 1 ^) }Z = ° (A10) 
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where the fields A®(xi) and ip^(x\) have been replaced by the derivatives of the generating functional with respect 
to the sources J af *(xi) and rj 7 (x\) and each of the subscripts a,j3 and 7 marks the components of color, flavor and 
spinor. Differentiating Eq. (A10) with respect to the source rj p (yi) and then setting 77 = 77 = ^ = ^ = 0, we obtain [5] 

[(id Xl - mi) a7 + (r>v)^ ^JL-j ]S F ( Xl - Vl ) J lp = S ap SHx 1 - yi )Z[J] (All) 

where 

Sf{xi - iSff^SnM ^=^=° (A12) 

is the quark propagator in the presence of source J. When the source J is set to be zero, noticing 

A;(x t I Xl , yi ) lp = -^—. S F {x 1 -y 1 ) J lp | J=0 (A13) 

here i = 1, 2, Eq. (All) straightforwardly goes over to the equation written in Eq. (2.10). 

Here, we think, it is appropriate to write down general definitions of some kinds of Green functions in the presence 
of source J. These Green functions will be encountered in later derivations and no longer defined repeatedly later on 
for saving pages. First, we write the remaining fermion propagators 

SAX2 ~ ^ " iS % (x 2 )5rf M (AU) 

C *, y7 ^ 2 Z[J,7y,77,g,^] 

and 

Myi " y2) - = w« '--^= o ■ (A16) 

On differentiating the above propagators with respect to the sources J atl ( Xi ), ... and J b "(yj), we get the following 
Green functions: 

K'.^.'X X //,.... I x uVl ) J ap - 777^7^ - 777^ (A17) 



^:±:{-n- I ^2,2/2)^ = , XTa " u ^,v, XT6 w., A -gg-(»2 - 2/2)^, (A18) 



and 



5 


S 


i8J a v( Xi )" 




5 


5 


iSJ a ^(x t )' 




5 


S 


i5J a ^{xi)' 


"i5J bv {y ) 


S 


5 


iSJ a ^(xi)' 


"iSJ b "( yj ) 



A; a ±;( Xi ,..., yj ,... I x u x 2 )i p = Tau , y .. . x Tbl/UiA ...S F ( Xl -x 2 ) J a0 (A19) 



A* a .'.',..'."(^, -,Vj,... I yuV2) J prT = -77^77rT-^77TT- 5 F(yi - Itt)^ (A20) 
where i, j = 1,2. The four-point Green function is defined, in the presence of source J, as 

G( Xl x 2 ;y u y 2 ) app „ - ^ a{xi)m w*MWM ^=^=° ' (A21) 

By successively differentiating the above Green function with respect to the sources J a>1 (xi), ... and J bu (yj),..., we 
obtain the Green functions like this 
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I •'•|,'-j://i-//ji;:,„. = /[U „ jt:r)! • • • 7-77^—- x 1 . z 2 : // : . ft • (A22) 

Taking differentiations of Eq. (A10) with respect to the sources rfp(x 2 ), ?7 P (j/i) and T}%{y 2 ) and noticing the equality 
in Eq. (A7) and the following nonvanishing derivatives 



it can be found that 



(id Xl - m 1 ) ai G(xix 2 ;yi,y2)j 0pa = 5 ap 5 i {x 1 - yi)S F (x 2 - y 2 )j a , 
+C aP S\ Xl -x 2 )S* F (yi -y2)% - {T^) ai Gl( Xl I x u x 2 ; yi ,y 2 )J 0pa [ > 

where we have set rj = rj = £ = £ = and the Green functions involved have been defined in Eqs. (A14), (A16), (A21) 
and (A22). When the source J is turned off, Eq. (A. 24) directly gives rise to the equation of motion in Eq. (2.24). 
A further differentiation of Eq. (A24) with respect to J b (x 2 ) gives 

{id Xl - mi) ai G b v {x 2 I xi.x 2 ;yi,y 2 ) J lfjpa = <5 ap (5 4 (a;i - yi )K cb (x 2 | x 2 ,y 2 )j a 
+C aP S 4 (x 1 -x 2 )At\x 2 I yi,y 2 ) J pa - {T a ^) ai G ab v {xi,x 2 \ xi,x 2 ;y u y 2 ) J lPprj 

where the Green functions have been defined in Eqs. (A18), (A20) and (A22). On letting the source J to vanish, Eq. 
(A25) will lead to the equation of motion in Eq. (2.30). 

Now, let us differentiate Eq. (A10) with respect to rf^(x 2 ). Once all the sources but J are set to be zero, we obtain 

[(id Xl - mi ) ai + (T a ^) ai J - }S* F (xix 2 )^ = {C- 1 ) a(i 5\x 1 - x 2 )Z[J] (A26) 

10J ^{Xl) 

where the propagator was defined in Eq. (A15). 

After differentiating Eq. (A26) with respect to J biy (x 2 ) and noticing the definition given in Eq. (A19), one gets 



(id Xl ~ m 1 ) ai A* l/ b (x 2 I xix 2 )Jp = (C 1 ) a pS i (x 1 - x 2 ) ^ SZ[J - 
-(T a ») aj A™ b (x 1 ,x 2 \x 1 ,x 2 ) J 1 p. 



'iSJ^W) (A27) 
When we set J = and considering i( ;jfi^ 2 ) |j=o= 0, Eq. (A27) simply gives the equation in Eq. (2.32). 



B. Equations of motion with respect to X2 

To derive the equations of motion with respect to x 2 for the Green functions , Let us take a derivative of the 
generating functional in Eq. (Al) with respect to the field i/jp(x 2 ). By the same procedure as described in Eqs. 
(A8)-(A10), it can be found that 

m x 2 ) + {i i8 X2 - m2hx + i^jJ^J^Z - (Bl) 

In the above derivation, the relation in Eq. (A7) has been used. Differentiating Eq. (Bl) with respect to r}%{y 2 ) and 
setting all the sources except for J to vanish, we have 

l(id X2 - m 2 ), x + (T a ")f>x iSJb t {x2) }S F (x 2 y 2 )l = 5 Pa 5\x 2 y 2 )Z[J] (B2) 

where S F (x 2 — y 2 )'{ a was defined in Eq. (A14). When we set J = and noticing 

K b (*2 I x 2 ,y 2 ) Xa = i§j I [x2) Sf^2 y 2 )L \j=o (B3) 

and the relation in Eq. (2.15), Eq. (B2) will become the equation in Eq. (2.14). 
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Upon differentiating Eq. (Bl) with respect to the sources rj a (xi), f] p {yi) and rj^(y 2 ) and setting all the sources but 
J to be zero subsequently, it can be found that 

(id X2 - m 2 ){3\G{x l x 2 \y ll y 2 )l Xp(y = Sp a S 4 (x 2 - y 2 )S F (x 1 - yi ) J ap 
+C al3 5 4 (x 1 -x 2 )S* F (yi -y<i) J pa - (T V )a\G h v {x 2 \ xi,x 2 ;yi,y 2 ) J aXpa 

where G h v {x 2 \ xi t x 2 ;yi,y 2 ) 

aXpa was defined in Eq. (A22). When the source J is turned off, Eq. (B4) just leads to 

the equation in Eq. (2.25). 



C. Equations of motion with respect to yi 

In order to find the equations with respect to y\ for the Green functions, it is necessary to differentiate the generating 
functional in Eq. (Al) with respect to the field ip p (yi)- Following the procedure stated in Eqs. (A8)-(A10), it is easy 
to get 

iW + [«K + - ^r Pj ^]}z = (ci) 

On differentiating the above equation further with respect to rj a (xi) and then keeping the source J to be nonvanishing 
only, we find 

S F ( Xl - yi)^ T [(*K + mi) Tp - l6J i [yi) ^)r P ] = -5 ap S\ Xl yi )Z[J] (C2) 

When the source J is turned off , noticing 

A£(m I x uyi ) aT = . g S F ( Xl -yi)L \ j= (C3) 

which is defined in Eq. (A17) and the relation 

A»(yi I x uVl ) aT (T a ^) Tp = J d i z 1 S F (x 1 - *i)ar£(*i,yi) Tp (C4) 

which is easy to be proved, Eq. (C2) just gives rise to the equation shown in Eq. (3.3). 

On differentiating Eq. (CI) with respect to rj a (xi), rj c p{x 2 ) and rj%{y 2 ) and then setting ^ = r/ = ^ = ^ = 0,we have 

G(x 1 x 2 ;y 1 ,y 2 )i fjTa (i < d yi + m x ) Tp = -5 ap 6 i {x 1 - yi)S c F (x 2 - y 2 ) J p<J 

-C pa S\y 1 -y 2 )S F ( Xl - X 2)^ + G-(y 1 \ Xl , X2] y 1 ,y 2 ) J al3Ta ^ afl )r P 1 ' 

where the propagators were defined in Eqs. (A14) and (A15), while the Green function G"(j/i | x\,x 2 ;y\,y 2 ) J was 
given in Eq. (A. 22). Clearly, the above equation directly leads to the equation in Eq. (3.5) when J is set to vanish. 
Taking the differentiation of Eq. (C5) with respect to J afl (xi), we get 

G°(xi I x 1 x 2 \y ll y 2 ) J af}Ta {id yi +m 1 ) Tp = -5 ap 5 4 {x 1 - yi )A c p a (xi \ x 2l y 2 )j ia , , 

-C prT S 4 (y 1 -y 2 )A; a (x 1 \ x u x 2 ) J a? + Gf v {xr, yi \ x 1 ,x 2 ;y 1 ,y 2 )i f3Ta (r b -)r P 1 ' 

where the Green functions on the RHS of the above equation have already defined in Eqs. (A18), (A19) and (A22). In 
the case of the source being absent, we get from Eq. (C6) the equation of motion in Eq. (4.4). While, differentiating 
Eq. (C6) with respect to J b "{x 2 ) further, one can get 

Gf l/ (x 1 ,x 2 I x l x 2 ;y l ,y 2 ) J a0Ta {idy 1 +m 1 ) Tp = -6 Q(9 <S 4 (a;i - yi)A°f(x 1 ,x 2 \ x 2 ,y 2 ) J f}a , , 

-C pa S 4 (y 1 -y 2 )A;f(x 1 ,x 2 \ x u x 2 ) J af3 + Gf v c K (x u x 2 ; Vl \ x u x 2 ; yi ,y 2 )ip Ta {T™) Tp 1 ' 

where all the Green functions were already defined in Eqs. (A18), (A19) and (A22). The equation given above 
straightforwardly gives rise to equation (4.22) when the source is set to be zero. 

Now we turn to derive an equation of motion obeyed by the propagator S F {y\ — y 2 ). By differentiating Eq. (CI) 
with respect to rj^(y 2 ) and then setting 77 = 77 = ^ = ^ = 0, it will be found 
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S* F (Vi ~ yiYrA^y, + m 1 ) Tp = C pa 8\yi - y 2 )Z[J} + At\ yi I y u y 2 )U^)r P (C8) 
where S* F (yi — y-i) J Ta and A* 6 (yi | y\,y<i) J Ta were denned in Eqs. (A16) and (A20). Considering 

b S 

A* (yi | yx,yi) r a = iSJ bv( y .y ?F(vi - yi) J Ta \j=o (C9) 

and the relation 

A* b (yi I yi, y 2 )r. (r bl/ ) rp = ^ d 4 z{s F ( Zl - y2)r^(z 1 ,y 1 ) Tp (cio) 

which is easy to be proved by using Eqs. (2.7) and (2.11), we see, Eq. (C8) isl immediately converted to thw equation 
inEq. (3.7)'. 

Taking a differentiation of Eq. (C8) with respect to J afl (x 2 ), we obtain 

A*°(x 2 I yi,y 2 )i a (u) yi +m 1 ) Tp = C pa 8 i {y 1 - y 2 ) a jfJ^ 2) 
+A* l 2 b (x 2 ,y 1 | y u y 2 )i a {T b ») Tp 

where the Green functions were defined in Eq. (A20). On setting J = 0, the above equation goes over to equation 
(4.39). 

D. Equations of motion with respect to y2 

For deriving the equations of motion with respect to y 2 for the Green functions, we need to differentiate the 
generating functional in Eq. (Al) with respect to ip^(y 2 ). By the procedure shown in Eq. (A8)-(A10), it can be found 



6 r/'*7T \ /T^> v \ $ 

■x TT~{ i( ld V2 + m 2)S* - r ) Sa —— nJ 



+ 71^7— M*dy. + ^2)6, - (T V ,, Tft ,^, , ]}Z = 0. (Dl) 



Then, we differentiate the above equation with respect to r]^(x 2 ) and subsequently set all the sources, except for J, 
to be zero. As a result , we get 

S F {x 2 - y 2 )j 5 [(d y2 + m 2 ) Sa l5J t {m) (r bi V] = -5 f 3„5 4 (x 2 y 2 )Z[J}. (D2) 
On setting J — and noticing 

At b (V2 | X 2 ,V2hs = i5 jt { y 2) S F^ - V*)U \J=0 ( D3 ) 

and 

A?(lte I x 2 , y2 )p S {T bu ) Sa = J d 4 z 2 S F (x 2 - z 2 ) (iS Y, c {z 2 - y 2 ) Su , (D4) 

we obtain the equation in Eq. (3.4) from Eq. (D2). 

To derive Eq. (4.8), we need to differentiate Eq. (D2) with respect to J aM (xi). The result is 

A™(xi | x 2l y 2 ) J fjx (d V2 + m 2 ) Xa = Sp a 5 4 (x 2 - y 2 ) i<5 jfi^ l} 
+A™ b ( Xl ;y 2 \x 2 ,y 2 y^(T b,/ ) Xa 



(D5) 



where the Green functions were represented in Eq. (A18). On setting J — 0, the above equation just leads to the 
equation in Eq. (4.8). 

Upon differentiating Eq. (D5) with respect to J b,/ (x 2 ), it is found that 



S 2 Z[J] 

j - -vpa" - y 2 /JsjaJl " 

+A c p a ^{x 1 ,x 2 ;y 2 \ x 2 , y 2 ) J ^ x (T CK ) x<T 



A p a J'(xi,x 2 | x 2 ,y 2 )j x (d y2 +m 2 ) Xa = -5 fia 5\x 2 - y 2 ) mKF ^± w ^ ) 
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where the Green's functions were also represented in Eq. (A18). In the case of J = 0, considering that 

5 2 Z[J] 



«^ (ai )^(* a -) |j=0=<A - (!Bl - a!a) (D?) 

which is the gluon propagator, the above equation straightforwardly gives rise to equation (4.26). 

Now we turn to derive the equation of motion for the four-point Green function. After completing the differentiations 
of Eq. (Dl) with respect to rj a {xi),rfp{x 2 ) and r] p (yi) and then letting r] = rj = £ = £ = 0, we have 



G(x 1 x 2 ;y 1 ,y 2 )^ /3p5 (id y2 + m 2 ) Sa = -5p a 5 {x 2 -y2)S F {xi -yi) 
-C pa d 4 (y 1 - y 2 )Sp(x 1 - x 2 ) J afj + G b u {y 2 \ x 1 ,x 2 ;y 1 ,y 2 )^ p5 (T v ) Sa 



(D8) 



where G b v {y 2 \ X\,x 2 ;yi,y 2 )^ g was defined in Eq. (A22). It is clear that when J = 0, the above equation leads to 
the equation written in Eq. (3.10). 

Going on further, it is seen that the differentiation of Eq. (D8) with respect to J a,1 (yi) gives 

I x 1 x 2 \y ll y 2 y arjp& {id y2 +m 2 ) s<y = -6p a 6 4 (x 2 - y 2 )A a p ( yi | a^I/i)^ ^ 
-C pa 5 4 { Vl -y 2 )A* p a (y 1 \ x u x 2 ) J afi + Gf v {yx, y 2 \ x 1 ,x 2 ;yi,y2)Zp pS (r v )& a 

where the Green functions were defined in Eqs. (A17), (A19) and (A22). Obviously, when the source J is turned off, 
the equation shown above goes over to equation (3.14). 

Differentiating Eq. (D8) with respect to J afl (xi), we obtain an equation like this 



G°(xi | xi x 2 ;yi,y 2 y a/Sp5 (idy 2 + m 2 ) S a = -5^5 4 (x 2 -y 2 )A%(xi \xi,yi)i 



■bu 



P 



-{C 1 ) pa S 4 (y 1 -y 2 )A* p a (x 1 \ x u x 2 ) J a o + G ab v {x x ; y 2 x u x 2 ; yi ,y 2 )Z 5 {T ) s , 



(D10) 



where the Green functions can also be found in Eqs. (A17), (A19) and (A22). Letting the source to be zero, the 
equation shown above will lead to the equation in Eq. (4.9). 

Taking a further differentiation of Eq. (D10) with respect to J bl/ (yi), it is seen that 

GfAxiWi I x 1 x 2 ;y 1 ,y 2 )j lflpS (i*d y2 +m 2 ) Sa = -6 f3tT S 4 (x 2 - y 2 )A ab v {x x ; y x \ x 1 ,y 1 )i p 
HC- 1 ) pa 5\ yi - y2 )Klf(x l] y l \x ll x 2 )i p + G^{x 1 -y 1 ,y 2 \ 

When the source is turned off, the equation above just gives the equation in Eq. (4.11). 
Now, let us take a differentiation of Eq. (D10) with respect to J bv (x 2 ). The result is 

G ab v {xi,x 2 | x l x 2 ;y l ,y 2 ) J al3pS {id y2 +m 2 ) 5(T = -Sp a S 4 (x 2 - y 2 )A ab v { Xl , x 2 \ x i,yi)t p , m2 s 
+ (C- 1 ) pa S 4 (y 1 -y 2 )A* p a J>(x 1 ,x 2 \ x u x 2 ) J af) + G a p b u c K (x u x 2 ; y 2 x u x 2l y u y 2 ) J aPpS {T CK ) Sl y 

where the Green functions can be read from Eqs. (A17), (A19) and (A22). When the source J is absent, the above 
equation will be converted to the equation in Eq. (4.28). A further differentiation of the above equation with respect 
to J CK {y\) leads to 

G P t c K (x ll x 2 ;y 1 \ x 1 x 2 ;yi,y 2 ) J af}pS {idy 2 +m 2 ) Sa - 
= SpJ 4 (x 2 -y 2 )A^ c K (xi,a;2;j/i | xi,yi)i p 

+(C- 1 ) P J 4 (yi-y 2 )A;t K c (x 1 ,x 2 ,y 1 \ x u x 2 ) J af3 ^ D16 > 
+G a p b u c ^(x 1 ,x 2 ;y 1 ,y 2 \ x u x 2 ; yi , y2)i f3pS (T C ) Sa 

where the Green functions can also be read off from Eqs. (A17), (A19) and (A22). When the source is turned off, 
Eq. (D13) will go over to the equation in Eq. (4.30). 

To derive the equations in Eqs. (3.11) and (3.17), it is necessary to differentiate Eq. (Dl) with respect to r) p {y\). 
On completing this differentiation, we find 

S* F (yi - y2) J pS [(dy 2 + m 2 )sa l§J t [v2) F bV )^ = C P j\yi V2)Z[J\ (D14) 

where we have set r] = rj = £ = £ = 0. Noticing 
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and 



b 5 

A* (y 2 | 2/1, y 2 ) P 5 = jjjw - I J=o ( D1 5) 



Al b (y 2 1 yi, y 2 )p 5 (r b ' y )^ = y ^2^(2/1 - *s)£s c (*2 - y 2 ) s „, (D16) 

equation (D14) will be brought to the equation in Eq. (3.11). 

By differentiating Eq. (D14) with respect to J aM (yi) and noticing the definition given in Eq. (A20), we obtain 



A* a (yi | yi,y2) J pS (d V2 + m 2 ) Sa = C p J 4 ( yi - y 2 ) tS jal J { l l} _ 

T* ab / I \ / /t$> v \ \ UL ' ) 

+ A M „ (2/1 , 2/2 I 2/1, Vi) p & ( r )Sa- 



Equation (3.17) directly follows from the above equation when we set J = 0. On taking the derivative of Eq. (D14) 
with respect to J afJ -(x 2 ) , it is clearly seen that 

A*"(x 2 I yi,y 2 ) J pS (dy 2 +m 2 )sa = C pa 5 4 ( yi - y 2 ) w jfifl a) 
+A M „ (ar 2 ;j/2 | 2/1,2/2)^^ ) 5ct 

where the Green functions were all shown in Eq. (A20). After setting J = 0, the equation above becomes the equation 
in Eq. (4.41). Differentiating the above equation further with respect to J b "{yi), we get 



A*f {x 2 - Vl I yi,y 2 y p5 (dy 2 +m 2 ) Sa = C pa 5 i {y 1 - 2/2) ^ J a M( t 2 ^ [ /] l ,, (yi) j 
[x 2 ;yi,y 2 \ yi,y 2 ) u pS (i ) Sa 



where the Green functions were represented in Eq. (A20). When the source is set to vanish and noticing the definition 
in Eq. (D7) for the gluon propagator, the above equation just gives the equation in Eq. (4.42). 
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